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ABSTRACT 

We analyze a 19-night photometric search for transiting extrasolar planets in the open cluster NGC 
1245. An automated transit search algorithm with quantitative selection criteria finds six transit 
candidates; none are bona fide planetary transits. We characterize the survey detection probability 
via Monte Carlo injection and recovery of realistic limb-darkened transits. We use this to derive upper 
limits on the fraction of cluster members with close-in Jupiter-radii, Rj, companions. The survey 
sample contains ~870 cluster members, and we calculate 95% confidence upper limits on the fraction 
of these stars with planets by assuming the planets have an even logarithmic distribution in semimajor 
axis over the Hot Jupiter (HJ - 3.0< P/day <9.0), Very Hot Jupiter (VHJ - 1.0< P/day <3.0), and 
an as of yet undetected Extremely Hot Jupiter (EHJ - PRoche < P <1.0 day) period ranges. For 
1.5 Rj companions we limit the fraction of cluster members with companions to <1.5%, <6.4%, and 
<52% for EHJ, VHJ, and HJ companions, respectively. For 1.0 Rj companions, we find <2.3% and 
<15% have EHJ and VHJ companions, respectively. We do not reach the sensitivity to place any 
meaningful constraints on 1.0 Rj HJ. From a careful analysis of the random and systematic errors 
of the calculation, we show that the derived upper limits contain a ±7%^° relative error. For similar 
photometric noise and weather properties as this survey, observing NGC 1245 twice as long results 
in a tighter constraint on HJ companions than observing an additional cluster of similar richness as 
NGC 1245 for the same length of time as this survey. If 1% of stars have 1.5 Rj HJ companions, we 
expect to detect one planet for every 5000 dwarf stars observed for a month. To reach a ^ 2% upper 
limit on the fraction of stars with 1.5 Rj companions in the 3.0< P/day <9.0 range, we conclude a 
total sample size of ~ 7400 dwarf stars observed for at least a month will be needed. Results for 1.0 
Rj companions, without substantial improvement in the photometric precision, will require a small 
factor larger sample size. 

Subject headings: open clusters and associations:individual (NGC 1245) — planetary systems — meth¬ 
ods: data analysis 


1. INTRODUCTION 

Extrasolar planet detections and analysis of the non¬ 
detections further our knowledge of the planet forma¬ 
tion process and contribute to an empirical determina¬ 
tion of the typical planetary system. These empirical 
constraints will eventually decide the ubiquity or rarity of 
planetary bodies in the Universe. A variety of techn iques 
exist to detect extrasolar planets ilPerrvmanll200(Tll . and 
there are currently 168 extrasolar planet discoveries'^. 

Most of the extrasolar planets have been discovered us¬ 
ing the radial velocity technique. Radial velocity detec¬ 
tions indicate that 1.2% ± 0.3% of FGK main-sequence 
stars in the solar neighborhood have a “ Hot” Jupiter- 
mass planet (HJ) orbiting within 0.1 AU llMarcv et al.l 
l2?io^ . At this small separation from the central star, 
the high temperatures and low disk c olumn densities pre¬ 
vent in situ formation of HJ planets llBodenheimer et al.l 
|2003). Several mechanisms exist to exchange angular 
momentum between the protoplanet and natal disk, en¬ 
abling the protoplanet to migrate from a more likely 
formation separation (several AU) to within 0.1 AU 
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llTerQuem et al.ll200r)t) . Due to tidal circularization, HJ 
have nearly circular orbits with an observed median ec¬ 
centricity, < e >= 0.07, whereas planets with larger sep¬ 
arations have a median eccentricity of < e >= 0.25. 

In addition to the detection statistics and planet prop¬ 
erties, the extrasolar planet detections indicate several 
physical relationships between the stellar host properties 
and the frequency of extrasolar planets. The most strik¬ 
ing of these is that the probability for hosting an extraso¬ 
lar planet increases rapidl y with stellar metal abun dance, 
consistent with P oc llFischer fc Valentill2005li . The 
frequency of planets may also depend o n the s tellar mass. 
iButler et al.l ^|2004^ and iBonfils et alJ l|2005ll point out 
that there exists a deficit of Mj planets orbiting M 
dwarf stars. However, the increasing number of short- 
period Neptune-mass planets being discovered around M 
dwarfs suggests that the overall frequency of planets (of 
all masses) orbiting M dwarfs may be similar to FGK 
dwarfs, but the typical planet mass is less, thereby es¬ 
caping detection given the det ection limitations of the 
current radial velocity surveys l|Bonfils et al.ll2?)f)^ . Ad¬ 
ditionally, none of the M dwa rfs harboring planets are 
metal rich IBonfils et all (120051) . 

A coherent theory of planet formation and survival re¬ 
quires not only reproducing the physical properties of the 
planets, but reproducing any trends in the physical prop¬ 
erties on the host environment. Despite the knowledge 
and constraints on extrasolar planets that radial velocity 
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surveys provide, radial velocity surveys have their limi¬ 
tations. The high resolution spectroscopic requirements 
of the radial velocity technique limit its use to the solar 
neighborhood and orbital periods equivalent to the life¬ 
time of the survey. A full consensus of the planetary for¬ 
mation process requires relying on additional techniques 
to detect extrasolar planets in a larger variety of condi¬ 
tions prevalent in the Universe. 

For instance, microlensing surveys are sensitive 
to extrasolar planets orbiting stars in the Galactic 
disk and bulge with distances of many kpc away 
ijMao fc PacvzhskilllQ^ iGould fc Loeblll992(). Two ob¬ 
jects consistent with Jupiter-mass companio ns have been 
detected via the microl ensing technique llBond et al.l 
l2?inllTTdalski et a,lj[2f)f)l . Additional information is ob¬ 
tained from studying the microlensing events that did not 
result in extrasolar planet detections. Microlensing sur¬ 
veys limit the fraction of M dwarfs in the Galactic bulge 
with Jupiter-mass com panions orbiting be t ween 1.5 AU 
to 4 AU to < 33% llAlbrow et alJ I2001t IGaudi et al.1 

l2?inl . 

Although limited to the solar neighborhood, attempts 
to directly image extrasolar planets are sensitive to plan¬ 
ets with semimajor axis beyond 20 AU. The light from 
the parent star limits detecting planets interior to the 
seeing disk. Adaptive optics observations of young (~ 1 
Myr) stars provide the best opportunity to directly im¬ 
age extrasolar planets since the young planets are still 
relatively bright while undergoing a rapid, cooling con¬ 
traction. Although the interpretation relies on theoreti¬ 
cal modeling of these complex planetary objects, three 
sources in nearby star forming regions have been de¬ 
tected whose broad-band colors and spectra are con¬ 
sistent wfth__those_ex£ected_fromJM2_Ju£iter-maas_^- 
jects llNeuhauser. et al.l l200,'Tt iGhauvin et al.l l2005albll . 
The contrast ratios necessary for extrasolar planet de¬ 
tection are difficult to reach, and results for detecting 
higher mass brown dwarfs are more complete. An anal¬ 
ysis of the Gornell High-Order Adaptive Optics Survey 
(GHAOS) derives a brown dwarf companion upper limit 
of 10% orbiting between 25 and 100 AU of the parent star 
llGarson et aljl2?)(i^ . iMcGarthv fc Zuckermanl ll2004ll es¬ 
timate 1% ± 1% of G,K, and M stars have brown dwarf 
companions orbiting between 75 and 300 AU, but this 
estimate may not account for the fu ll range of orbital in - 
clination and eccentricities possible llGarson et al.ll2fiC)^ . 
At greater separations, > 1000 AU, brown dwarf com¬ 
panions to F-MO main-sequenc e stars appear to be as 
common as stellar companions llGizis et al ] |200!I1 . 

After the radial velocity technique, the transit tech¬ 
nique h as had the most suc cess in detecting extrasolar 
planets llKonacki et al.ll2005ll . The transit technique can 
detect Rj transits in any stellar environment where <1% 
photometry is possible. Thus, it provides the possibility 
of detecting extrasolar planets in the full range of stellar 
conditions present in the Galaxy: the Solar neighbor¬ 
hood, the thin and thick disk, open clusters, the halo, 
the bulge, and globular clusters are all potential targets 
for transit surveys. A major advantage of the transit 
technique is the current large-format mosaic CCD im¬ 
agers which provide multiplexed photometric measure¬ 
ments with sufficient accuracy across the entire field of 
view. 

The first extrasolar planet detections via the tran¬ 


sit technique began with the candidate list p rovided by 
the OGLE collaboration l|Udalski et al.ll20d^ . However, 
confirmation of the transiting extrasolar planet candi¬ 
dates requires radial velocity observations. Due to the 
well known equation-of-state competition between elec¬ 
tron degeneracy and ionic Coulomb pressure, the radius 
of an object becomes insensitive to mass across the en- 
tire range from below Mj to the hydrogen-burning limit 
ijChabrier fc Baraffel200nll . Thus, objects revealing a Rj 
companion via transits may actually have a brown-dwarf 
mass companion when followed up with radial velocities. 
This degeneracy is best illustrated by the planet-sized 
brow n dwarf companion to OGLE-TR-122 (|Pont et alJ 
ioo^. The first radia l-velocity confirmations of planets 
disco vered by transits llKonacki et alJ 200 ^lBouchv et alJ 
l 2 ?iril provided a first glimpse at a population of massive, 
very close-in planets with P < 3 days and Mp > Mj 
(“Very Hot Jupiters” - V HJ) that had not b een seen 
by radial velocity surveys. IGaudi et al.l ll2005ll demon¬ 
strated that, after accounting for the strong sensitivity of 
the transit surveys to the period of the planets, the tran¬ 
sit detections were likely consistent with the results from 
the radial velocity surveys, implying that VHJs were 
intrinsically very rare. Subs equently, in a metall icity- 
biased radial velocity survey, iBouchv et al.1 ll2005bll dis¬ 
covered a VHJ with P = 2.2 day around the bright star 
HD189733 that also has observable transits. 

Despite the dependence of transit detections on radial 
velocity confirmation, radial velocity detections alone 
only result in a lower limit on the planetary mass, and 
thus do not give a complete picture of planet formation. 
The mass, radius information directly constrains the the¬ 
oretical models, whereas either parameter alone does lit¬ 
tle to further constrain the impor tant physical processes 
that shape the planet properties llGuillotll2005ll . For ex¬ 
ample, the mass-radius relation for extrasolar planets 
can constrain the si ze of the rocky core present (e.g., 
iLauEfhlin et alJl2005ll . Also, the planet transiting across 
the face of its parent star provides the exciting poten¬ 
tial to probe the planetary atmosph eric absorption lines 
against the stellar spectral features llCharb onneaii et alJ 

l2002tlDeming et alJ2005atlNarita et al.l2005ll . Or. in the 

opposite case, emission from the planetary atmosphere 
can be detected when the planet orbits behind the par - 
ent star llCharbonneaii et alJ2005l:fDeming et al.l2005bll . 

Despite these exciting results, the transit technique is 
significantly hindered by the restricted geometrical align¬ 
ment necessary for a transit to occur. As a result, a 
transit survey necessarily contains at least an order of 
magnitude more non-detections than detections. In ad¬ 
dition, null results themselves can provide important con¬ 
straints. For example, the null result in the globular clus¬ 
ter 47 Tucanae adds important empirical constraints to 
the trend of increasing probability of hav ing a planetary 
companion with increas ing metallicity llGi11i1a,nd et alJ 
1200(1 l^ntos et ani2004ll . Thus, understanding the sen- 
sitiyity of a giyen transit suryey, i.e. the expected rate of 
detections and non-detections, takes on increased impor¬ 
tance. Seyeral studies have taken steps toward sophis¬ 
ticated Monte Carlo calculations to quantify detection 
probabilities in a transit suryey |lGil lilandet_aJjl 2 QQd; 

2005HHood et aLlboOfiH . Unfortunately, these studies do 

not fully characterize the sources of error and systemat- 
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ics present in their analysis, and therefore the reliability 
of their conclusions is unknown. Furthermore, essentially 
all of the previous studies have either (1) not accurately 
determined the number of dwarf main-sequence stars in 
their sample, or (2) made simplifying assumptions which 
may lead to misestimated detection probabilities, or (3) 
contained serious conceptual errors in the procedure with 
which they have determined detection probabilities, or 
(4) some combination of the above. 

As a specific and important example, most studies do 
not apply identical selection criteria when searching for 
transits amongst the observed light curves and when re¬ 
covering injected transits as part of determining the sur¬ 
vey sensitivity. Removal of false-positive transit candi¬ 
dates arising from systematic errors in the light curve 
has typically involved subjective visual inspections, and 
these subjective criteria have not been applied to the re¬ 
covery of injected transits when determining the survey 
sensitivity. This is statistically incorrect, and can in prin¬ 
ciple lead to overestimating the survey sensitivity. Even 
if identical selection criteria are applied to the original 
transit search and in determining the survey sensitivity, 
some surveys do not apply conservative enough selections 
to fully eliminate false-positive transit detections. 

In this paper, we address these shortcomings of pre¬ 
vious studies in our analysis of a 19-night photometric 
search for transiting extrasolar planets in the open clus¬ 
ter NGC 1245. An automated transit search algorithm 
with quantitative selection criteria finds six transit can¬ 
didates; none are bona fide planetary transits. We de¬ 
scribe our Monte Carlo calculation to robustly determine 
the sensitivity of our survey, and use this to derive upper 
limits on the fraction of cluster members with close-in, 
Jupiter-radii, Rj, companions. 

Leading up to the process of calculating the upper 
limit, we develop several new analysis techniques. First, 
we develop a differential photometry method that au¬ 
tomatically selects comparison stars to reduce the sys¬ 
tematic errors that can mimic a transit signal. In ad¬ 
dition, we formulate quantitative transit selection crite¬ 
ria, which completely eliminate false positives due to sys¬ 
tematic light-curve variability without human interven¬ 
tion. We characterize the survey detection probability 
via Monte Carlo injection and boxcar recovery of tran¬ 
sits. Distributing the Monte Carlo calculation to mul¬ 
tiple processors enables rapid calculation of the transit 
detection probability for a large number of stars. 

The techniques developed here enable combining re¬ 
sults from transit surveys in a statistically meaningful 
way. This work is part of the Survey for Transiting 
Extrasolar Planets in Stellar Systems (STEPSS). The 
project concentrates on stellar clusters since they pro¬ 
vide a large sampl e of stars of homogeneo us metallicity, 
age, and distance lIBurke et al.ll2003l 1200411 . Overall, the 
project’s goal is to assess the frequency of close-in extra¬ 
solar planets around main-sequence stars in several open 
clusters. By concentrating on main-sequence stars in 
open clusters of known (and varied) age, metallicity, and 
stellar density, we will gain insight into how these var¬ 
ious properties affect planet formation, migration, and 
survival. 

The survey characteristics and the photometric proce¬ 
dure are given in in We explain the automated algo¬ 
rithm to calculate the differential light curves and de¬ 


scribe the light curve noise properties in In 0we de¬ 
scribe our impl ementation of the bo x-fitting least squares 
(B LS) m ethod llKovacs et alJl200'^ for transit detection. 
In Sin we present a thorough discussion of the quan¬ 
titative selection criteria for transit detection, followed 
by a discussion of the objects with sources of astrophys- 
ical variability that meet the selection criteria in We 
outline the Monte Carlo calculation for determining the 
detection probability of the survey in m We present up¬ 
per limits for a variety of companion radii and orbital 
periods in in A discussion of the random and system¬ 
atic errors present in the technique is given in 3H1 We 
compare the final results of this study to our expected 
detection rate before the survey began and discuss the 
observations necessary to reach sensitivities similar to 
radial velocity detection rates in m Finally, ^1 ni briefly 
summarizes this work. 

2. OBSERVATIONS AND DATA REDUCTION 
2 .1. Observations 

We observed NGC 1245 for 19 nights between 24 Oct. 
and 11 Nov. of 2001 using the MDM 8K mosaic imager 
on the MDM 2.4m Hiltner telescope. The MDM 8K im¬ 
ager consists of a 4x2 array o f thinned, 2048x4096, SITe 
ST002A CCDs e ro 113 12001 . This instrumental setup 
yields a 26'x26' field of view and 0.36" per pixel reso¬ 
lution in 2x2 pixel binning mode. Tabled has an entry 
for each night of observations that shows the number of 
exposures obtained in the cousins J-band filter, median 
full-width-at-half-maximum (FWHM) seeing in arcsec- 
onds, and a brief comment on the observing conditions. 
In total, 936 images produced usable photometry with a 
typical exposure time of 300 s. 

2 .2. Data Reduction 

We use the IRAF® CCDPROC task for all CCD pro¬ 
cessing. The read noise measured in zero-second images 
taken consecutively is consistent with read noise mea¬ 
sured in zero-second images spread through the entire 
observing run. Thus, the stability of the zero-second 
image over the course of the 19 nights allows median 
combining 95 images to determine a master, zero-second 
calibration image. For master flat fields, we median com¬ 
bine 66 twilight sky flats taken throughout the observing 
run. We quantify the errors in the master flat field by 
examining the night-to-night variability between individ¬ 
ual flat fields. The small-scale, pixel-to-pixel variations 
in the master flat fields are ~ 1%, and the large-scale, 
illumination-pattern variations reach the 3% level. The 
large illumination-pattern error results from a sensitiv¬ 
ity in the illumination pattern to telescope focus. How¬ 
ever, such large-scale variations do not affect differential 
photometry with proper reference-star selection (as de¬ 
scribed in 0. 

To obtain raw instrumental photometric measure¬ 
ments, we employ an automated reductio n pipeline that 
uses the DoPHOT PSF fitting package llSchecter et alJ 
rr^ . Comparable quality light curves resulted from 
photometry via the DAQ PHQT/ALLFRAME, PSF- 
htting photometry packages (|StetsoiJll987HStetson et all 

^ IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 













4 


Burke et al. 


Il99i^ in the background limited regime. DoPhot per¬ 
forms slightly better in terms of rms scatter in the differ¬ 
ential light curve in the source-noise limited regime. The 
photometric pipeline originated from a need to produce 
real-time photometry of microlensing events in order to 
search for anomalies indicatin g the presence of an extra¬ 
solar planet around the lens llAlbrow et al linn). This 
study uses a variant of the original pipeline developed at 
The Ohio State University and c urrently in use b y the 
Microlensing Follow Up Network llYoo et al.ll2004^ . 

In brief, the pipeline takes as input a high signal-to- 
noise (S/N) “template” image. A first pass through 
DoPHOT identifies the brightest, non-saturated stars on 
all the images. Using these bright-star lists, an auto¬ 
mated routine (J. Menzies, priyate communication) de¬ 
termines the geometric transformation between the tem¬ 
plate image and all the other other images. A second, 
deeper pass with DoPhot on the template image identifies 
all the stars on the template image for photometric mea¬ 
surement. The photometric procedure consists of trans¬ 
forming the deep-pass star list from the template image 
to each frame. These transformed positions do not yary 
during the photometric solution. Next, an automated 
routine (J. Menzies, priyate communication) determines 
an approximate value for the FWHM and sky as required 
by DoPHOT. Finally, DoPHOT iteratively determines a 
best-fit, 7-parameter analytic PSF and uses this best-fit 
PSF to determine whether an object is consistent with a 
single star, double star, galaxy, or artifact in addition to 
the photometric measurement of the object. 

3. DIFFERENTIAL PHOTOMETRY 

In its simplest form, differential photometry involves 
the use of a single comparison star in order to remove the 
time variable atmospheric e xtinction signal from the raw 
photometric measurements llKieldsen fc Frandsenlll993l . 
The process of selecting comparison stars typically con¬ 
sists of identifying an ensemble of bright, isolated stars 
that demonstrate long term stability over the course of 
the observations llGilliland fc Brownlll988ll . This proce¬ 
dure is sufficient for studying many variable astrophys- 
ical sources where several percent accuracy is typically 
adequate. However, after applying this procedure on a 
subset of the data, systematic residuals remained in the 
data that were similar enough in shape, time scale, and 
depth to the expected signal from a transiting compan¬ 
ion to result in large number of highly-significant false 
positive detections. 

Removing ;^0.01 mag systematic errors resembling a 
transit signal requires a time consuming and iterative 
procedure for selecting the comparison ensemble. Addi¬ 
tionally, a comparison ensemble that successfully elim¬ 
inates systematic errors in the light curve for a partic¬ 
ular star fails to eliminate the systematic errors in the 
light curve of a different star. Testing indicates each 
star has a small number of stars or even a single star to 
employ as the comparison in order to reduce the level 
of systematics in the light curve. On the other hand, 
Poisson errors in the comparison ensemble improve as 
the size of the comparison ensemble increases. Addi¬ 
tionally, the volume of photometric data necessitates an 
automated procedure for deciding on the “best” possible 
comparison ensemble. Given its sensitivity to both sys¬ 
tematic and Gaussian noise and its efficient computation. 


we choose to minimize the standard deviation around the 
mean light curve level as the figure of merit in determin¬ 
ing the “best” comparison ensemble. 

3.1. Differential Photometry Procedure 

We balance improving systematic and Poisson errors in 
the light curve using the standard deviation as the figure 
of merit by the following procedure. The first step in de¬ 
termining the light curve for a star is to generate a large 
set of trial light curves using single comparison stars. 
We do not limit the potential comparison stars to the 
brightest or nearby stars, but calculate a light curve us¬ 
ing all stars on the image as a potential comparison star. 
All comparison stars have measured photometry on at 
least 80% of the total number of images. A sorted list of 
the standard deviation around the mean light-curve level 
identifies the stars with the best potential for inclusion 
in the comparison ensemble. Calculation of the stan¬ 
dard deviation of a light curve involves 3 iterations elim¬ 
inating 3-standard-deviation outliers between iterations. 
However, the eliminated measurements not included in 
calculation of the standard deviation remain in the final 
light curve. 

Beginning with the comparison star that resulted in 
the smallest standard deviation we continue to add in 
comparison stars with increasingly larger standard de¬ 
viations. At each epoch, we median combine the results 
from all the comparison stars making up the ensemble af¬ 
ter removing the average magnitude difference between 
target and comparison. We progressively increase the 
number of stars in the comparison ensemble to a maxi¬ 
mum of 30, calculating the standard deviation of the light 
curve between each increase in the size of the compari¬ 
son ensemble. The final light curve is determined using 
the comparison ensemble size that minimizes the stan¬ 
dard deviation. Less than 1% of the stars result in the 
maximum of 30 comparison stars. The median number 
of comparison stars is 4, with a modal value of 1. The 
distribution of comparison stars has a standard deviation 
around the median of 4. The fact that the standard devi¬ 
ations of the majority of stars is minimized using a single 
comparison star emphasizes the importance of consider¬ 
ing all stars as possible comparisons in order to minimize 
systematic errors and achieve the highest possible accu¬ 
racy. 

3.2. Comparison to a Similar Algorithm 

Independent of this study, iKovacs et al.l ll200f)ll devel¬ 
oped a generalized algorithm for eliminating systematic 
errors in light curves that shares several basic properties 
with the method we have just presented. They agree 
with the conclusion that optimal selection of comparison 
stars can eliminate systematics in the light curve. They 
also use the standard deviation of the light curve as their 
figure of merit (see their Equation 2). However, their 
more general implementation allows for the comparison 
star to have a real-valued, linear correlation coefficient 
(ci in their Equation 1) in the differential photometry, 
whereas the implementation outlined here forces binary 
values, 0 or 1, for the linear correlation coefficient. They 
solve for the real-valued, linear correlation coefficients by 
minimization of the standard deviation via matrix alge¬ 
bra, whereas the method given here relies on brute force 
minimization of the standard deviation. 
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A thorough comparison of the performance between 
these methods has not been done. However, we empha¬ 
size that our algorithm found that the modal number of 
stars in the comparison ensemble is a single star. Their 
algorithm restricts the comparison ensemble to a subset 
of the available stars. The restricted comparison ensem¬ 
ble may not capture the systematics present in a light 
curve. However, their real-valued, linear correlation coef¬ 
ficients may provide the degree of freedom lacking in the 
algorithm of this study necessary to cancel the systemat¬ 
ics. Both algorithms possess an important caveat. The 
figure of merit cannot distinguish between improvements 
in the Poisson error or systematic error and therefore 
does not guarantee optimal elimination of the system¬ 
atic deviations. 

3.3. Additional Light-curve Corrections 

Although our procedure for optimally choosing com¬ 
parison stars succeeds in dramatically reducing system¬ 
atics in the light curves, we find that some additional sys¬ 
tematic effects nevertheless remain. We introduce several 
additional corrections to the light curves to attempt to 
further reduce these effects. 

In good seeing, brighter stars display saturation ef¬ 
fects. Whereas, in the worst seeing, some stars display 
light-curve deviations that correlate with the seeing. To 
correct for these effects, we fit a two-piece, third-order 
polynomial to the correlation of magnitude versus see¬ 
ing. The median seeing separates the two pieces of the 
fit. We first fit the good-seeing piece with the values of 
the polynomial coefficients unconstrained. We then fit 
the poor-seeing piece, but constrain the constant term 
such that the fit is continuous at the median seeing. How¬ 
ever, we do not constrain the first or higher order deriva¬ 
tives to be continuous. In performing this fit, we excise 
measurements from the light curve that would lead to 
a seeing-correlation correction larger than the standard 
deviation of the light curve. We use this two-piece fit to 
correct the measurements. 

Measurements nearby bad columns on the detector also 
display systematic errors that are not removed by the 
differential photometry algorithm. Thus, measurements 
when the stellar center is within 6 pixels of a bad column 
on the detector are eliminated from the light curve. 

The final correction of the data consists of discarding 
measurements that deviate by more than 0.5 mag from 
the average light curve level. This prevents detection of 
companions with radii > 3.5i?j around the lowest mass 
stars of the sample. 

3.4. Light-curve Noise Properties 

Figure^shows the logarithm of the standard deviation 
of the light curves as a function of the apparent /-band 
magnitude. Calculation of the standard deviation in¬ 
cludes one iteration with 3-standard-deviation clipping. 
To maintain consistent S/N at fixed apparent magnitude, 
the transformation between instrumental magnitude to 
apparent /-band magnitude only includes a zero-point 
value, since including a color term in the transformation 
results in stars of varying spectral shape and thus varying 
S/N in the instrumental I band having the same appar¬ 
ent /-band magnitude. Each individual CCD in the 8K 
mosaic has its own zero point, and the transformation is 
accurate to 0.05 mag. 



Fig. 1. — Shows the logarithm of the light-curve standard devia¬ 
tion as a function of the apparent 7-band magnitude (points). The 
depths of transits due to a 1.0 and 1.5 Rj companion assuming 
the star is a cluster member are shown as dashed lines. The solid 
lines show photometric noise models that match the empirically 
determined noise properties. 


One CCD has significantly better noise properties than 
the others as evidenced by the second sequence of points 
with improved standard deviation at fixed magnitude. 
The instrument contains a previously unidentified prob¬ 
lem with images taken in the binning mode. The data 
was taken with 2x2 native pixels of the CCD array binned 
to one pixel on readout. During readout, the control sys¬ 
tem apparently did not record all counts in each of the 
four native pixels. However, the single CCD with im¬ 
proved noise properties does not suffer from this problem 
whereas all the other CCDs do. Subsequent observations 
with large positional shifts allow photometric measure¬ 
ments of the same set of stars on the affected detectors 
and unaffected detector. Performing these observations 
in the unbinned and binning modes confirms that on the 
affected detectors, 50% of the signal went unrecorded by 
the data system. This effectively reduces the quantum 
efficiency by half during the binned mode of operation 
for seven of the eight detectors. 

The two solid lines outlining the locus of points in Fig¬ 
ure ^ provide further evidence for the reduction in quan¬ 
tum efficiency. These lines represent the expected noise 
due to a a source-noise limited error, a term that scales as 
a background-noise limited error, and 0.0015 mag noise 
floor. We determine the lower line by varying the area 
of the seeing disk and the flat noise level until the noise 
model visually matches the locus of points for the de¬ 
tector with the lower noise properties. Then, the upper 
line results from assuming half the quantum efficiency 
of the lower noise model while keeping the noise floor 
the same. The excellent agreement between the higher 
noise model and the noise properties of the remaining 
detectors strongly supports the conclusion that half of 
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the native pixels are not recorded during readout. This 
readout error could introduce significant errors in the 
limit of excellent seeing. However, only 4% of the pho¬ 
tometric measurements have FWHM<2.5 binned pixels. 
Thus, even in the binning mode, we maintain sufficient 
sampling of the PSF to avoid issues resulting from the 
readout error. 

The different noise properties between detectors does 
not complicate the analysis. The transit detection 
method involves merit criteria (see that nat¬ 

urally handle data with varying noise properties. Other 
than reducing the overall effectiveness of the survey, the 
different noise properties between the detectors does not 
adversely affect the results in any way. 

In addition to the empirically determined noise prop¬ 
erties, DoPhot returns error estimates that on average 
result in reduced = 0.93 for a flat light-curve model. 
The average reduced x^ for all the detector agree within 
10%. Scaling errors to enforce reduced x^ = 1-0 for each 
detector independently has a negligible impact on the 
results, thus we choose not to do so. 

The upper and lower dash lines in Figure ^ show the 
transit depth assuming the star is a cluster member for 
1.5 and 1.0 Rj companions, respectively. In Figure ^ 
3671 stars have light curves with a standard deviation 
less than the signal of a transiting Rj companion. 

4. TRANSIT DETECTION 

In the previous section, we describe a procedure for 
generating light curves that reduces systematic errors 
that lead to false-positive transit detections. However, 
systematics nevertheless remain that result in highly sig¬ 
nificant false-positive transit detections. This section de¬ 
scribes the algorithm for detecting transits and meth¬ 
ods for eliminating false positives based on the detected 
transit properties. There are two types of false-positives 
we wish to eliminate. The first is false-positive transit 
detections that result from systematic errors imprinted 
during the signal recording and measurement process. 
The second type of false-positive results from true astro- 
physical variability that does not mimic a transit signal. 
For example, sinusoidal variability can result in highly 
significant detections in transit search algorithms. We 
specifically design the selection criteria to trigger on tran¬ 
sit photometric variability that affects a minority of the 
measurements and that are systematically faint. How¬ 
ever, the selection criteria do not eliminate false-positive 
transit signals due to true astrophysical variability that 
mimic the extrasolar planet transit signal we seek (graz¬ 
ing eclipsing binaries, diluted eclipsing binaries, etc.). 

For detecting transits w e employ the box-fitt ing least 
squares (BLS) method of iKovacs et al.l l|2002ll . Given 
a trial period, phase of transit, and transit length, the 
BLS method provides an analytic solution for the tran¬ 
sit depth. We show in the appendix the equivalence of 
the BLS method to a x^ minimiz ation. Instead of us ing 
the Signal Residue (Equation 5 in iKovac s et al.ll2 002ll or 
Signa l Detection Efficiency fEniiation 6 in lKo^cT^ abl 
[M for quantifying the significance of the detection, 
we use the resulting improvement in x^ of the solution 
relative to a constant flux fit, as outlined in the appendix. 

This section begins with a discussion of the parameters 
affecting the BLS transit detection algorithm. We set 
the BLS algorithm parameters by balancing the needs 


of detecting transits accurately and of completing the 
search efficiently. The next step involves developing a 
set of selection criteria that automatically and robustly 
determines whether the best-fit transit parameters result 
from bona fide astrophysical variability that resembles a 
transit signal. A set of automated selection criteria that 
only pass bona fide variability is a critical component 
of analyzing the null-result transit survey and has been 
ignored in previous analyses. 

Due to the systematic errors present in the light curve, 
statistical significance of a transit with a Gaussian noise 
basis is not applicable. In addition, the statistical signif¬ 
icance is difficult to calculate given the large number of 
trial phases, periods, and inclinations searched for tran¬ 
sits. Given these limitations, we empirically determine 
the selection criteria on the actual light curves. Although 
it is impossible to assign a formal false alarm probability 
to our selection criteria, the exact values for the selection 
criteria are not important as long as the cuts eliminate 
the false positives while still maintaining the ability to 
detect Rj objects, and identical criteria are employed in 
the Monte Garlo detection probability calculation. 

4.1. BLS Transit Detection Parameters 

The BLS algorithm has two parameters that determine 
the resolution of the transit search. The first parameter 
determines the resolution of the trial orbital periods. The 
BLS algorithm (as implemented bv IKovacs et al ] l200l 
employs a period resolution with even frequency inter¬ 
vals, ~ where Pi is the previous trial orbital 

period, P 2 is the subsequent (longer) trial orbital period, 
and r] determines the frequency spacing between trial 
orbital periods. During implementation of the BLS al¬ 
gorithm, we adopt an even logarithmic period resolution 
by fractionally increasing the p eriod, P 2 = Pi x (1 + rf). 
The original implementation bv IKovacs et al.l ll2002li for 
the orbital-period spacing is a more appropriate proce¬ 
dure, since even frequency intervals maintain constant or¬ 
bital phase shifts of a measurement between subsequent 
trial orbital periods. The even logarithmic period reso¬ 
lution we employ results in coarser orbital phase shifts 
between subsequent trial orbital periods for the shortest 
periods and increasingly finer orbital phase shifts toward 
longer trial orbital periods. Either period-sampling pro¬ 
cedure remains valid with sufficient resolution. We adopt 
77 = 0.0025, which, given the observational baseline of 19 
days, provides <10% orbital phase shifts for orbital pe¬ 
riods as short as 0.5 day. 

The second parameter of the BLS algorithm deter¬ 
mines the resolution in orbital phase by binning the 
phase-folded light curve. Binning of the data in orbital 
phase drastically improves the numerical efficiency, but 
not wi thout loss in determ ining the correct transit prop¬ 
erties. IKovacs et al ](|2002) give a thorough examination 
of how the sensitivity in recovering transits varies with 
orbital-phase binning resolution. To search for transit 
candidates in the light curves we adopt A^bins = 400 
orbital-phase bins. We verify with tests that the above 
parameters accurately recover boxcar signals in the light 
curves. After injection of boxcar signals in the light 
curves, we calculate the x^ of th® solution returned by 
the BLS method with the x^ of the injected model. Tests 
show that the BLS method with the above parameters 
return a x^ within 30%, and typically much better, of 
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the injected model’s ^. 

4.2. Selection Criteria 

We apply the BLS method following the description in 
the previous section to search for transit candidates in all 
6787 stars with light curves. A visual inspection of a light 
curve folded at the best-fit transit period can generally 
be used to discriminate between bona fide astrophysi- 
cal variability and a false positive arising from system¬ 
atic errors. However, a proper statistical assessment of 
the sensitivity of a transit search requires that the exact 
same set of selection criteria that are applied to cull false 
positives are also applied when assessing the detection 
probability via, e.g., Monte Carlo injection and recovery 
of artificial signals. Due to the large number of artificial 
signals that must be injected to calculate the detection 
probability properly, using selection criteria based on vi¬ 
sual inspection of light curves is practically very difficult 
or impossible. Therefore, quantitative, automated detec¬ 
tion criteria that mimic the visual criteria must be used. 

We employ four selection criteria that eliminate all 
false detections while still maintaining the ability to de¬ 
tect Rj companions. These four selection criteria con¬ 
stitute cuts on (1) the improvement in of a transit 
model over a constant flux model, (2) the ratio between 
the of the best-fit transit model and the of the 
best-fit anti-transit model, (3) the fraction of Ay^ from 
a single night, and (4) the transit period. 

The first of the selection criteria is a cut on Ax^, the 
improvement in between a constant flux model and a 
transit model. The Ay^ is s imilar to the Signal Residue, 
SR, of iKovacs et al.l ll2002ll : we derive Ay^ and its re¬ 
lation to SR in the appendix. We prefer Ax^ over SR 
as the former allows a direct comparison of the transit 
detection significance between light curves with different 
noise properties. Given the correlated systematics in the 
data, we cannot rely on analytical formulations with a 
Gaussian statistics basis for the statistical significance of 
a particular Ay^ value. We empirically determine a cut 
on Ay^ in combination with the other selection criteria 
to fully eliminate false detections. For a transit detection 
we require Ay^ > 95.0. As shown in the appendix, this 
selection criterion corresponds to a S/N~10 transit de¬ 
tection. Figure 121 shows the Ay^ of the best fit transit for 
all light curves along the x-axis. The vertical line desig¬ 
nates the selection criteria on this parameter. Even with 
such a strict threshold, there are still a large number of 
false positives that pass the Ay^ cut. 

Systematic variations in the light curves that are char¬ 
acterized by small reductions in the apparent flux of star 
that are coherent over the typical time scales of plane¬ 
tary transits can give rise to false-positive transit detec¬ 
tions. However, under the reasonable expectation that 
systematics do not have a strong tendency to produce 
dimming versus brightening of the apparent flux of the 
stars, one would expect systematics to also result in a 
false-positive ‘anti-transit’ (brightening) detections. Fur¬ 
thermore, most intrinsic variables can be approximately 
characterized by sinusoids, which will also result in sig¬ 
nificant transit and anti-transit detections. On the other 
hand, a light curve with a true transit signal and insignif¬ 
icant systematics should produce only a strong transit 
detection, and not a strong anti-transit detection. 

Thus, the ratio of the significance of the best-fit tran- 
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Fig. 2.— The small points show as a function of for 
the resulting best-fit transit parameters in all light curves. Here 
Ay^ and Ay^ are the improvement between the fiat light-curve 
model and the best-fit transit and anti-transit model, respectively. 
The dotted vertical line shows the Ay^ = 95.0 selection bound¬ 
ary. The solid diagonal line shows the Ay^/Ay^ = 2.75 selection 
boundary. Objects in the lower right corner pass both selection 
criteria. The green diamonds show values of Ay^ and Ay?, for the 
six transit candidates. The blue stars show the recovered values 
of Ay^ and Ay^ for the four light curves with injected transits 
shown in Figure^ The label next to the blue stars corresponds to 
the label in the upper right corner of each panel in Figure|7| These 
curves were created by injecting transits into the same light curve. 
The blue star labeled 0 shows the values of Ay^ and Ay_ for this 
light curve before the example transits were injected. 


sit signal relative to that of the best-fit anti-transit sig¬ 
nal provides a rough estimate of the degree to which 
a detection has the expected properties of a bona fide 
transit, rather than the properties of systematics or si¬ 
nusoidal variability. In other words, a highly significant 
transit signal should have a negligible anti-transit sig¬ 
nal, and therefore we require the best-fit transit to have 
a greater significance than the best-fit anti-transit. We 
accomplish this by requiring transit detections to have 
Ay^/Ay^ > 2.75, where Ay?, is the y^ improvement of 
the best-fit anti-transit. For a given trial period, phase 
of transit, and length of transit, the BLS algorithm re¬ 
turns the best-fit transit without restriction on the sign 
of the transit depth. Thus, the BLS algorithm simulta¬ 
neously searches for the best-fit transit and anti-transit, 
and so determining Ay?, has no impact on the numerical 
efficiency. 

Figure 121 shows the Ay?, of the best fit anti-transit ver¬ 
sus the Ay^ of the best-fit transit for our light curves. 
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The diagonal line demonstrates the selection on the ra¬ 
tio Ax^/Ax?. = 2.75. Objects toward the lower right 
corner of this Figure pass the selection criteria. The ob¬ 
jects with large Ax^ typically have correspondingly large 
Ax-■ This occurs for sinusoidal-like variability or strong 
systematics that generally have both times of bright and 
faint measurements with respect to the mean light-curve 
level. 

Requiring observations of the transit signal on sepa¬ 
rate nights also aids in eliminating false-positive detec¬ 
tions. We quantify the fraction of a transit that occurs 
during each night based on the fraction of the transit’s 
X^ signihcance that occurs during each night. The pa¬ 
rameters of the transit allow identihcation of the data 
points that occur during the transit. We sum the indi¬ 
vidual Xi = values for data points occurring 

during the transit to derive XtotJ where rui is the light 
curve measurement and Oi is its error. Then we calcu¬ 
late the same sum for each night individually. We denote 
this Xfeth night- We identify the night for which xLhnight 
contributes the greatest fraction of Xtot; we call this 
fraction / = xlth night/xL- Finally, we require / < 0.65. 
This corresponds to roughly seeing the transit one and a 
half times assuming all observations have similar noise. 
Alternatively, this criterion is also met by observing 2/3 
of a transit on one night and 1/3 of the transit on a 
separate night, or observing a full transit on one night 
and 1/6 of transit on a separate night with three times 
improvement in the photometric error. Figure |21 shows 
/ versus the best-fit period for all the light curves. The 
horizontal line designates the selection on this parameter. 

The red points in Figure El show objects that pass the 
Ax^ > 95.0 selection. We find that most are clustered 
around a 1.0 day orbital period. A histogram of the best- 
fit transit periods amongst all light curves reveals a high 
frequency for 1.0 day and 0.5 day periods. Visual inspec¬ 
tion of the phased light curves reveals a high propensity 
for systematic deviations to occur on the Earth’s rota¬ 
tional period and 0.5 day alias. We do not fully under¬ 
stand the origin of this effect, but we can easily conjec¬ 
ture on several effects that may arise over the course of 
an evening as the telescope tracks from horizon to hori¬ 
zon following the Earth’s diurnal motion. In order to 
eliminate these false positives, we apply as our fourth 
selection criteria a cut on the period. Specihcally, we 
require transit detections to have periods that are not 
within 1.0 ±0.1 and 0.5 ±0.025 day. The horizontal lines 
designate these ranges of discarded periods. 

5. TRANSIT CANDIDATES 

Six out of 6787 stars pass all four selection criteria. 
All of these stars are likely real astrophysical variables 
whose variability resembles that of planetary transit light 
curves. However, we find that none are bona fide plan¬ 
etary transits in NGC 1245. After describing the prop¬ 
erties of these objects we will describe the procedure for 
ruling out their planetary nature. Figure El shows the 
phased light curves for these six stars. Each light-curve 
panel in Figure 0 has a different magnitude scale with 
fainter flux levels being more negative. The upper left 
corner of each panel gives the detected transit period as 
given by the BLS method. The upper right corner of each 
panel gives an internal identification number. The panels 
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Fig. 3.— The black points show / as a function of the best-fit 
transit orbital period, where / is the fraction of the total 
provement with the best-fit transit model that comes from a single 
night. The objects that pass the Ax^ > 95.0 selection criteria are 
shown as red points. The horizontal line shows the / = 0.65 se¬ 
lection boundary. The vertical lines denote orbital period regions 
avoided due to false-positive transit detections. The blue stars and 
green diamonds are the same as in Figure 1^ 


from top to bottom have decreasing values in the ratio 
between the improvement of a transit and anti-transit 
model, Ax^/Ax?_. 

TableEllists the properties and selection criteria values 
for the stars shown in Figure El The green diamonds in 
Figures El and El represent the selection criteria for the 
six transit candidates. T he photometric an d positional 
data in TableElcome from iBurke et al.l ll2004ll . The Xmem 
entry in Table El measures the photometric distance of a 
star from the isochrone that best fits the cluster CMD. A 
lower value of this parameter means a star has a position 
in the CMD closer to the main sequence. Heavy points 
in FigureEldenote stars with Xmem < 0.04, and we desig¬ 
nate these stars as potential cluster members. Based on 
XmeniJ Star 20513 and star 70178 have photometry consis¬ 
tent with cluster membership, thus we also l ist the phys¬ 
ical p arameters of those stars in Table El iBurke et alJ 
^I2004^ details the procedure for determining the physi¬ 
cal parameters of a star based solely on the broad-band 
photometry and the best-fit cluster isochrone. However, 
the validity of the stellar physical parameters only ap¬ 
plies if the star is a bona fide cluster member. 

Figure El shows a finding chart for each star with a 
light curve in Figure El The label in each panel gives the 
identification number, and the cross indicates the corre- 
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Fig. 4.— The points show the change in magnitude as a func¬ 
tion of orbital phase for all stars that meet the transit candidate 
selection criteria. Negative values for A mag are toward fainter 
flux levels. The phased period is given in the upper left corner of 
each panel, and the number in the upper right corner of each panel 
gives the internal identiflcation number. 



Fig. 5.— The CMD of the cluster NGC 1245. Potential clus¬ 
ter members having Xmem ^ 0.04 are shown with heavy points). 
Objects that exceed the selection criteria for transit detection are 
given as open diamonds. 


spending object. Star 20274 is not centered in the finding 
chart because it is located near the detector edge. The 
field of view of each panel is 54". North is toward the 
right, and East is toward the bottom. The panels for 
stars 20065, 20398, and 20513 (located near the cluster 



Fig. 6.— Finding charts for transit candidates with light curves 
shown in Figure HI Each panel is 54" on a side. North is toward 
the right, and East is toward the bottom. 


center) provide a visual impression of the heaviest stel¬ 
lar crowding encountered in the data. Figure 0 shows 
the V and B — V CMD of the cluster field as given in 
iBurke et al.l ll2004ll . The large open stars denote the lo¬ 
cations of the objects that exceed the transit selection 
criteria. 

5.1. Consistency of Transit Parameters with Cluster 
Membership 

Only stars 20513 and 70718 have values consis¬ 
tent with cluster membership. Additionally, the transit 
depth in both stars indicates potential for having a Rj 
companion. However, qualitatively, in each case the tran¬ 
sit duration relative to the orbital period is too long to be 
a true planetary companion to a cluster main-sequence 
star. We can use our knowledge of the physical properties 
of the parent stars to quantitatively rule out planetary 
companions. We do this by comparing an estimate of the 
stellar radius derived from the CMD to an independent 
estimate of a lower limit on the stellar radius derived 
from the properties of the light curve. In both cases we 
find that the stellar radii derived from the CMD are well 
below the lower limit on the stellar radius based on the 
light curve. 

To derive a lower limit on the stellar radius 
from the light curve, we b uild on the work by 
iSeager fc Mallen-Ornelaj ll2003H . They provide a purely 
geometric relationship between the orbital semimajor 
axis, a, and stellar radius, i?*, for a light curve with 
a given period, P, depth of transit, AF, and total dura¬ 
tion of the transit (first to fourth contact), t, assuming a 
circular orbit (see their Equation 8). By assuming a cen¬ 
tral transit (impact parameter 6 = 0), we transform their 
equality into a lower limit. Using Kepler’s Third Law, 
assuming that the mass of the companion is much smaller 
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than the mass of the star, and assuming the duration of 
the transit is much smaller than the period (r <C P), we 
find, 


P 1/3(1 + v^)’ 


( 1 ) 


where P* is in AU, M* is in units of Mq and r and P 
are in years. 

Parameters on the right hand side of the above equa¬ 
tion contain substantial uncertainties. Replacing the 
parameters by their maximum plausible deviation from 
their measured values in such a manner as to decrease P* 
increases the robustness of the lower limit. The orbital 
period determination has the largest uncertainty. Tests 
of recovering transits in the light curves reveal a 10% 
chance for the BLS method to return an orbital period, 
P', at the 1/2P and 2P aliases of the injected orbital 
period, and a < 1% chance of detecting the 1/3P and 
3P aliases. Misidentification of the correct orbital pe¬ 
riod results from gaps in the observing window function. 
Replacing P in the above equation with 3P', where P' is 
the orbital period returned by the BLS algorithm, pro¬ 
vides the maximum plausible deviation of this quantity 
and increases the robustness of the lower limit. In addi¬ 
tion, the stellar mass determination based on the CMD 
potentially has contamination from a binary companion. 
Thus, we replace M* with 0.5M', where M'^ is the stel¬ 
lar mass estimate from the CMD. We do not modify r 
or AP. For the cases considered here, AP ^ 1, and the 
term 1 -|- VAP ~ 1 in Equation ^ Therefore the pre¬ 
cise value of AP has little effect on the resulting limit on 
P*. The BLS algorithm fits a boxcar transit model to 
the light curve via a minimization. Since, in the limit 
of zero noise, any non-zero boxcar height fit to a tran¬ 
sit can only result in an increasing when the length 
of the boxcar exceeds the length of the transit, r un¬ 
derestimates the true transit length. Making the above 
replacements the lower limit on the stellar radius is. 


(M;/MQ)i/3(^/lday) 
(PVlday)i/3(H-v/^) 


( 2 ) 


For star 20513, the above equation requires P* > 
1.04Pq if the star is a cluster member. Fits to the CMD 
yield a stellar radius P* = O.SOPq. The lower limit for 
star 70718 is P* > 0.82Pq, whereas the CMD yields 
P* = 0.56Pq. Clearly both stars lack consistency be¬ 
tween the stellar radius based on the CMD location and 
the stellar radius based on the transit properties. 

The transiting companions to 20513 and 70718 are 
also unli kely to be planets if the host stars are field 
dwarfs. iTinglev fc SacketH ll2005ll provide a diagnos¬ 
tic to verify the planetary nature of a transit when 
only the light curve is a vailable. The diagnostic rip of 
ITinglev fc SacketH ll200,5ll compares the length of the ob¬ 
served transit to an estimate of the transit length derived 
by assuming a main-sequence mass-radius relation for the 
central star. By assuming a radius of the companion of 
Rp = l.OPj, we find r]p = 4.0 and 3.8 for 20513 and 
70718, respectively. Values otr]p <1 correspond to plan¬ 
etary transits. Therefore, 20513 and 70718 are unlikely 
to host planetary companions with Rp < Rj if they are 
main-sequence stars. 

We note that our final a posteriori criterion with which 
we reject cluster transit candidates, namely the consis¬ 


tency between the radius of the parent star as estimated 
from the CMD and the radius as estimated from the light 
curve, is a conceptually different kind of selection crite¬ 
rion than those that we applied to all the light curves 
to arrive at our six transit candidates. The original four 
selection criteria were designed to detect bona fide astro- 
physical variability that resembles the signals from tran¬ 
siting planets, but does not necessarily arise from a tran¬ 
siting planetary companion. In principle, we could have 
included the radius consistency cut as an additional selec¬ 
tion criterion applied to all light curves. The motivation 
to do this would be that imposing this additional crite¬ 
rion might automatically remove some systematic false 
positives and so allow us to improve our efficiency by 
making the other selection criteria less stringent. We 
have found using limited tests that this is not the case. 
We therefore chose to leave the radius check as an a pos¬ 
teriori cut on the transit candidates. Nevertheless, ob¬ 
serving a cluster does provide an advantage over observ¬ 
ing field stars, as the additional constraint on the stellar 
radius from the cluster CMD provides a more reliable 
confir mation of the planetary n ature than the light curve 
alone ijTinglev fc Sackettll200,5(l . and furthermore allows 
a more accurate assessment of the detection probability. 

It is important to emphasize that all of the injected 
transits with which we compute the detection probability 
(§0 automatically pass the radius consistency criterion. 
A fraction of these will be recovered at periods that differ 
enough from the input period that by using the recovered 
period they will no longer satisfy the radius constraint. 
However, we find that this fraction is negligibly small. 

5.2. Individual Cases 

This section briefly discusses each object that met the 
selection criteria as a transit candidate but does not be¬ 
long to the cluster. The V-shaped transit detected in star 
30207 rules out a Rj companion. Transiting Rj compan¬ 
ions result in a flat bottomed eclipse as the stellar disk 
fully encompasses the planetary disk. A closer inspec¬ 
tion of the light curve also reveals ellipsoidal variations 
outside of the transit. This light curve matches the prop¬ 
erties of a grazing eclips e, which is a typical contaminant 
in transit searches le.g.. iBouchv et al]l200,5ll . 

The remaining stars have depths too large for a Rj 
companion and show evidence for secondary eclipses. Re¬ 
call that we eliminated data points with | Ato| >0.5 mag 
in the light curves. This eliminates the eclipse bottom 
for star 20065. Keeping all the data for star 20065 clearly 
reveals the characteristics of a detached eclipsing binary. 
The period BLS derives for star 20065 aligns the primary 
and secondary eclipses, and thus BLS-reported period is 
not the true orbital period. 

The eclipses in stars 20398 and 20274 do not perfectly 
phase up. This is because the resolution in period we 
used for the search prevents perfect alignment of the 
eclipses for such short periods. This effect is inconse¬ 
quential for detecting transiting planets as they all have 
orbital periods longer than 0.3 day. 

Finally, we note that other variables exist in the 
dataset. They were not selected because they do not 
meet the Ax^;j^/Ax^;jj_ selection criterion. A future 
paper will present variables that exist in this dataset 
using selection criteria more appropriate for identifying 
quasi-sinusoidal periodic variability (J. Pepper et ah, in 
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preparation). 

6. DETECTION PROBABILITY CALCULATION 

We did not detect any transit signals consistent with a 
Rj companion. To interpret this null result in terms of 
the frequency of planetary companions to stars in NGC 
1245, we develop a Monte Carlo detection probability 
calculation for quantifying the sensitivity of the survey 
for detecting extrasolar planet transits. The calculation 
provides the probability of detecting a transit in the sur¬ 
vey as a function of the companion semimajor axis and 
radius. In addition to the photometric noise and ob¬ 
serving window, the observed properties of the transit 
signal depend sensitively on the host mass, radius, limb- 
darkening parameters, and orbital inclination with re¬ 
spect to the line of sight. Without accurate knowledge 
of the stellar parameters, a detailed detection probabil¬ 
ity is not possible. This precludes analyzing stars not 
belonging to the cluster. Given the degeneracy between 
broad-band colors of dwarfs, subgiants, and giants, the 
stellar radius for most field objects cannot be determined 
from the CMD alone. Assuming all stars of a given color 
are dwarfs drastically overestim ates the number of actua l 
dwarf stars in a transit survey lIGould fc Morganll200,‘lll . 
The minimal expenditure of observational resources nec¬ 
essary for determining the stellar parameters for a cluster 
transit survey provides a significant advantage over tran¬ 
sit surveys of the field. 

Each star in the survey has a unique set of physical 
properties and photometric noise, thus we calculate the 
detection probability for all stars in the survey. This is 
the first study of its kind to do so. Given the detection 
probability for each star, the distribution of extrasolar 
planet semimajor axis, and frequency of extrasolar planet 
occurrence, the survey should have detected, 

Ndet = f*Y.Pdet,U (3) 

extrasolar planets, where the sum is over all stars in the 
survey, 

Rp^Pmem,idRpd(lj 

( 4 ) 

Rp is the extrasolar planet radius, a is the semimajor 
axis, /* is the fraction of stars with planets distributed 

according to the joint probability distribution of Rp and 
,2 

a, The Monte Garlo detection probability calcu¬ 

lation provides (a, i?p),the probability of detecting a 
transit in a given light curve. The term Pr,i(a, Rp) gives 
the probability for the planet to cross the limb of the host 
along the line of sight, and Pmem.i gives the probability 
the star is a cluster member. This framework for calcu¬ 
lating the e xpected detections of the survey follows from 
the work of iGaudi et al.1 ll2002ll . In the following subsec¬ 
tions we describe the procedure for calculating each of 
these probability terms. 

6.1. Calculating Rp) 

P^^i{a,Rp) is the probability of detecting a transit 
around the fth star of the survey averaged over the or¬ 
bital phase and orbital inclination for a given companion 


radius and semimajor axis. We begin this section with a 
description of the procedure for injecting limb-darkened 
transits into light curves for recovery. After injecting the 
transit, we attempt to recover the transit employing the 
same BLS algorithm and selection criteria as employed 
during the transit search on the original data. It is crit¬ 
ical to employ identical selection criteria during the re¬ 
covery as the original transit search since only then can 
we trust the robustness and statistical significance of the 
detection. The fraction of transits recovered for fixed 
semimajor axis and Rp determines P^. Next, we charac¬ 
terize the sources of error present in P^ and how we en¬ 
sure a specified level of accuracy. Finally, in this section 
we discuss the parallelization of the calculation to obtain 
Pe for all stars in the survey in a reasonable amount of 
time. 

In the appendix, we discuss th e importance of inject¬ 
ing realistic transits for recovery. iMandel fc Agoll ll2002fl 
provide analytic formulas for calculating realistic limb- 
darkened transits. We employ the functional form of a 
transit for a quadratic limb-dark ening law as given in 
Section 4 of iMandel fc Agoll ll200^. The quadr atic limb- 
darkening coefficients come from iGla.retl <1200(111 . Specifi¬ 
cally, we use the /-band limb-darkening coefficients using 
the ATLAS calculation for logg = 4.5, log[M/H]=0.0, 
and Vturb = 2 kms“^. 

We assume circular orbits for the companions. All 
known extrasolar planets to date that orbit within 0.1 
AU have eccentricities <0.3, and the average eccentric¬ 
ity for these planets is < e >= 0.07®. 

After injecting the transit, we employ the BLS algo¬ 
rithm to recover the injected transit signal using the 
selection criteria described in S21 For numerical effi¬ 
ciency, we relax the resolution of the BLS search param¬ 
eters. We adopt a fractional period step, ry = 0.004, and 
phase space binning, A^bins = 300. Despite the reduced 
resolution, higher resolution, converged solutions reveal 
only a 0.003 lower probability resulting from the adopted 
parameters. We correct all probabilities for this system¬ 
atic even though it is at an insignificant level compared 
to the other uncertainties. 

Figure|7|visualizes the injected transits with increasing 
degrees of significance from top to bottom. This Figure 
shows light curves with an injected transit phased at the 
period as returned from the BLS algorithm. The solid 
line illustrates the injected limb-darkened transit signal. 
The top two panels and the bottom two panels illustrate 
1.0 and 1.5 Rj companions, respectively. The transit 
recovery in the top panel barely meets the selection cri¬ 
teria, thus giving a visual impression for the sensitivity 
of the survey. The resulting selection criteria values after 
recovery of the injected transits are shown in Figures [3 
and 13 by the blue stars, and the labels next to the stars 
correspond to the panel label given in the upper right 
hand corner. The modeled transits shown in this Fig¬ 
ure are injected into the same light curve of a potential 
cluster member with V = 16.6 and rms scatter (before 
transit injection) of cr = 0.003. The blue star, labeled 0 
in Figures 13 and 13 represents the values for the selection 
criteria for found for this light curve before injecting the 
transits. 

As opposed to previous work, we carefully examine, 

® http://www.obspm.fr/encycl/catalog.html 
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Fig. 7.— Phased light curves showing the recovery of transits 
injected in the light curve by the Monte Carlo calculation. The 
injected limb-darkened transit signal is given by the solid line. The 
top two panels and bottom two panels show results for 1.0 and 1.5 
Rj companions, respectively. The transit recovery in the top panel 
barely meets the selection criteria and gives an impression for the 
sensitivity of the survey. The labels in the upper right corner of the 
panels correspond to the markers for the selection criteria values 
shown in Figures!^ and|3 


quantify, and control the uncertainties present in the cal¬ 
culation. During injection of a transit at fixed semima¬ 
jor axis, the transit can occur during any phase of the 
orbit. We use the following procedure to ensure that 
we inject enough trial transits at random orbital phases 
to yield convergence of Based on binomial statis¬ 
tics, the error in the resulting probability at fixed orbital 
period depends on the actual probability and the num¬ 
ber of trial transit phases, ^- Cact), 

where iVtriai is the number of trial transit phases and 
tact is the actual probability (unknown a priori). Main¬ 
taining the same error in the detection probability for 
differing tact requires a variable number of trial phases. 
For each semimajor axis, we first obtain an initial esti¬ 
mate for the probability, Cest, using fVtriai = 100. We 
then increase Vtriai until the probability converges to 
CTe = \/Vtriaieest(l “ Eest) < 0.02. The above proce¬ 
dure systematically overestimates Cact when tact ^ 0.95 
and systematically underestimates Cact when Cact ^ 0.05. 
However, these errors are of order the adopted = 0.02 
accuracy, and so we neglect them here. 

In addition to a random orbital phase, assuming a 
random orientation of the orbit requires taking into ac¬ 
count an even distribution in cos*, where i is the incli¬ 
nation of the orbit. Only a narrow range of inclinations, 
cos* < (i?* -I- Rp)/a, results in a transit. Thus, we in¬ 
ject the transit with an even distribution in cos* between 
0 < cos* < (i?* -I- Rp)/a. 

The previous discussion pertains to ensuring a pre¬ 
scribed accuracy at fixed semimajor axis. However, the 
expected detection rate also requires an integral over 
semimajor axis, which must be sampled at high enough 
resolution to ensure convergence of the integral. We 
calculate the probability at even logarithmic intervals, 


5 log a = 0.011 AU. In comparison to high-resolution, 
converged calculations, this semimajor axis resolution re¬ 
sults in an absolute error in the detection probability 
integrated over the semimajor axis of = 0.003. We 
inject transits with semimajor axis from the larger of 
0.0035 AU and 1.5i?* to 0.83 AU. The best-fit isochrone 
to the cluster CMD determines the parent star radius. 

Generating the light curve from the raw photometric 
measurements is numerically time consuming. Thus, we 
inject the transit after generating the light curve. This 
procedure has the potential to systematically reduce or 
even eliminate the transit signal, because generating the 
light curve and applying a seeing decorrelation tend to 
“flatten” a light curve. To quantify the significance of 
this effect, we inject transits in the raw photometric mea¬ 
surements before the light curve generation procedure on 
several stars in the sample. Comparing the detection 
probability obtained by injecting transits before light 
curve generation to the detection probability obtained 
by injecting the transit after light curve generation re¬ 
veals that injecting the transit after generating the light 
curve overestimates the detection probability by ^ 0.03. 
We decrease the calculated probability at fixed period by 
0.03 to account for this systematic effect. 

The 0.03 systematic overestimate in the detection 
probability becomes increasingly important for correctly 
characterizing the detection probability at long orbital 
periods. For instance, the detection probability for a star 
of median brightness will be overestimated by >15% for 
orbital periods > 4.0 day and 1.5 Rj companions if this 
systematic effect is not taken into account. The detection 
probability is overestimated by >50% for orbital periods 
>8.0 day without correction. The results for 1.0 Rj com¬ 
panions are even more severe. The detection probability 
would be overestimated by 50% for periods beyond 1.8 
day for a star of median brightness 

Based on the CMD of NGC 1245 
this study contains light curves for ~ 2700 stars consis¬ 
tent with cluster membership. Initially, we calculate the 
detection probability for 2 possible companion radii: 1.0 
and 1.5 Rj. For each star, on average we inject 50000 
transits for a single companion radius at 150 different 
semimajor axes. In total, we inject and attempt to re¬ 
cover ~ 2.7X 10® transits. Current processors allow injec¬ 
tion and attempted recovery on order of 1 s per transit. 
A single processor requires ~3000 days for the entire cal¬ 
culation. Fortunately, the complete independence of a 
transit injection and recovery trial allows parallelization 
of the calculation. We accomplish a parallel calculation 
via a server and client architecture. A server injects a 
transit in the current light curve and sends it to a client 
for recovery. 

Based on the computing resources available, we em¬ 
ploy two different methods for communication between 
the server and clients. Using a TCP/IP UNIX socket 
implementation for communication between the server 
and clients allows access to ~40 single-processor personal 
workstations connected via a local area network within 
the department of astronomy at The Ohio State Uni¬ 
versity. Additionally, the department of astronomy at 
The Ohio State University has exclusive access to a 48 
processor Beowulf cluster via the Cluster Ohio program 
run by the Ohio Supercomputer Center. The Message 
Passing Interface (MPI) libraries provide communication 


without correction. 
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Fig. 8.— Detection probability as a function of the orbital period 
is shown as the heavy solid line. This is a product of the proba¬ 
bility for a transit to occur (dash line) and the probability that 
an injected transit meets the selection criteria (light solid line). 
The panels from top to bottom show representative stars in order 
of increasing apparent magnitude. The left panels give results for 
a 1.5 Rj companion. The right panels give results for a 1.0 Rj 
companion. 


between the server and clients on the Beowulf cluster. 
A Beowoulf cluster belonging to the Korean Astronomy 
Observatory also provided computing resources for this 
calculation. C programming source code for either client- 
server communication implementation is available upon 
request from the author. 

The light solid line in Figure |H| shows the detection 
probability, P^{a,Rp), for three representative stars in 
order of increasing apparent magnitude from top to bot¬ 
tom and for the 2 companion radii, 1.0 and 1.5 Rj, on 
the left and right, respectively. In general, the proba¬ 
bility nears 100% completion for orbital periods <1.0 
day and then has a power law fall off toward longer or¬ 
bital periods. The falloff in the detection probability 
toward longer orbital periods partially results from the 
requirement of observing more than one transit. The 
large drop in the detection probability around 0.5 and 
1.0 day orbital periods results from the selection criteria 
we impose. The narrow, non-zero spikes in the detection 
probability near the 0.5 and 1.0 day orbital periods re¬ 
sult from injecting a transit at this period, but the BLS 
method returns a best-fit period typically at the ^0.66 
day alias. 

FigurelHlshows the detection probability with 3.3 times 
higher resolution in orbital period and a lower, 1%, error 
in the detection probability at fixed orbital period than 
the actual calculation. Thus, the figure resolves variabil¬ 
ity in the detection probability as a function of orbital 
period for probabilities >1%. However, such fine details 
have negligible impact on the results. 

6.2. Calculating Px^i{a, Rp) 

The probability for a transit to occur is Pt = {R* + 
Rp)/a. This transit probability assumes the transit is 
equally detectable for the entire possible range of orbital 


inclinations that geometrically result in a transit. As 
cosi for the orbit approaches (i?* -I- Rp)/a the transit 
length and depth decreases, degrading the transit S/N. 
We address this when computing P^ by injecting the tran¬ 
sit with an even distribution in cosi between the geo¬ 
metric limits for a transit to occur. Thus, Pt represents 
the overall probability for a transit with high enough in¬ 
clination to begin imparting a transit signal, while the 
detailed variation of the light curve signal for varying in¬ 
clination takes place when calculating P^. Pt is shown 
as the dashed light line in Figure |H1 The heavy solid line 
in Figure IHl is the product of Pe and Pt- 

6.3. Calculating Pmem 

The Monte Carlo calculation requires knowledge of the 
stellar properties, and the given properties are only valid 
if the star is in fact a bona fide cluster member. An 
estimate of the field-star contamination from the CMD 
provides only a statistical estimate of the cluster mem¬ 
bership probability. Based on th e study of the mass func- 
tion and field contamination in iBurke et alJ ll2004fl . we 
estimate the cluster membership probability, Pmem) as a 
function of stellar mass. In brief, we start with a subsam¬ 
ple of stars based on their proximity to the best-fit cluster 
isochrone (selection on Xmem < 0.04, see (El). This sam¬ 
ple contains V* ~ 2700 potential cluster members, and 
the heavy points in Figure [S] mark this cluster sample in 
the CMD. The best-fit isochrone allows an estimate of 
the stellar mass for each member of the cluster sample, 
and we separate the sample into mass bins. Repeating 
this procedure on the outskirts of the observed field of 
view, scaled for the relative areas, provides an estimate of 
the field-star contamination in a given mass bin. We fit 
Pmem, given in discrete mass bins, with a smooth spline 
fit for interpolation. 

The solid line in Figure |51 shows Pmem as a function 
of stellar mass. The corresponding probability is given 
on the right side ordinate. The clear histogram shows 
the distribution of the potential cluster members as a 
function of mass. The lower shaded histogram shows 
the product of the potential cluster members histogram 
and Pmem- This results in effectively ~870 clus¬ 

ter members in total. For reference, the corresponding 
apparent /-band magnitude is given along the top. 

7. RESULTS 

7.1. Results Assuming a Power-law Orbital-period 
Distribution 

The previous section describes the procedure for cal¬ 
culating the sensitivity of the survey to detect planetary 
companions as a function of semimajor axis. The re¬ 
sults from this calculation enable us to place an upper 
limit on the fraction of cluster members harboring close- 
in companions given the null result. However, calculating 
the upper limit over a range of orbital periods necessi¬ 
tates assuming a distribution of orbital periods for the 
planetary companions. Radial velocity surveys charac¬ 
terize the distribution of extrasolar planets in period as 
dn oc P~'^dP, with 0.7 < 7 < 1.0, corresponding to 
dn (X a~^da, with 0.5 < 0 < 1.0 ijSteDinski fc Blacfl 
l200ll iTabachnik fc Tremainel ^2002^ . These studies fit 
the entire range of orbital periods ranging from several 
days to several years. More recently, after an increase in 
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Fig. 9.— Distribution of the potential cluster members as a 
function of stellar mass {open histogram). The solid line shows the 
membership probability (right hand ordinate) as a function of stel¬ 
lar mass. The shaded histogram shows the product of the potential 
cluster member histogram and the cluster membership probability. 
The corresponding apparent 7-band magnitude is given along the 
top. 


the nu mber of extrasolar planet discoveries, lUdrv et al.l 
(l20(B confirm a shortage of planets with 10 < P < 100 
day orbits. Thus, the period distribution may take on 
different values of 7 in the P ^ 10 day and P > 100 day 
regimes. 

The initial extrasolar planet discoveries via the transit 
techn ique had periods less than 3.0 days l|Kona,cki et a,].l 
l 2 ?iol . The detection of these “Very Hot Jupiters,” 
contrasted with the results from radial-velocity surveys, 
which demonstrated a clear paucity of planets with P < 
3.0 days. After accounting for the strong decrease in 
sensi tivity of field trans it surveys with increasing pe¬ 
riod, iGaudi etlill l|2005j) demonstrated the consistency 
between the apparent lack of VHJ companions in the 
radial velocity surveys and their discovery in transit sur¬ 
veys. They further demonstrated that VHJ appear to be 
intrinsically much rarer than HJ (with 3 < P/day < 9. 
We will therefore treat VHJ as HJ distinct populations. 

Due to the incomplete knowledge of the actual period 
distribution of extrasolar planets and its possible depen¬ 
dence on the properties of the parent star, we provide 
upper limits assuming an even logarithmic distribution 
of semimajor axis. Thus, we assume a form of the joint 
probability distribution of the semimajor axis and Rp 
given by 

^ (5) 

where k is the normalization constant, S is the Dirac delta 
function, and Rp is the planet radius. We initially give 



Fig. 10.— Probability for transit detection as a function of the 
apparent 7-band magnitude assuming an even logarithmic distri¬ 
bution in semimajor axis from 1.0< P <3.0 day. The top panels 
assume Pmem = 1-0. The left panels show results for a 1.5 Rj 
companion. The right panels show results for a 1.0 Rj companion. 
The bottom panels are the same as the top panels, but they take 
into account the membership probability Anem- 


results for i?(, = 1.0 and 1.5 Rj. We follow IcTaudi et all 
and show results for HJ (3.0<P<9.0 day) and 
VHJ (1.0<P<3.0 day) ranges. In addition, we show re¬ 
sults for a more extreme population of companions with 
^Roche < P < 1.0 day, where PRoche is the orbital pe¬ 
riod at the Roche separation limit, which we designate 
as Extremely Hot Jupiter (EHJ). Assuming a negligible 
companion mass, the Roche period depends solely on the 
density of the companion. Jupiter, Uranus, and Neptune 
have nearly the same PRoche ~ 0.16 day. 

Figure M shows the probability for detecting a VHJ 
(1.0 day < P < 3.0 day) companion with an even log¬ 
arithmic distribution in semimajor axis as a function of 
apparent /-band magnitude. The left and right panels 
show results for a 1.5 and 1.0 Rj companion, respectively. 
The top panels of Figure EH show the probability for de¬ 
tecting an extrasolar planet, Pdet, assuming Pmem = 1 - 0 . 
The bottom panels show P^et after taking into account 
Pmem- The results for 1.0 Rj companions broadly scatter 
across the full range of detection probability. However, 
the 1.5 Rj companion results delineate a tight sequence 
in detection probability as a function of apparent mag¬ 
nitude. 

The 1.5 Rj companion signal lies many times above 
the rms scatter in the light curve (see Figure 0. Thus, 
a single measurement contributes a large fraction of the 
S/N required for detection. In this limit, the observing 
window function mainly determines the detection prob¬ 
ability, and as we show in H9.2I the result is similar to 
results obtained by the t heore tical detection probabil¬ 
ity framework of iGaii dl (| 2003 ). However, the 1.0 Rj 
companion transit signal comes closer to the detection 
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Fig. 11.— Upper limit (95% Confidence) on the fraction of stars 
in the cluster with companions for several ranges in orbital period 
assuming an even logarithmic distribution in semimajor axis. The 
solid lines show results for a 1.5 Rj companion. The dash lines 
show results for a 1.0 Rj companion. 


threshold. iPeoper fc Gaiidil ll2005ll describe the sensitiv¬ 
ity of a transit survey as a function of planet radius. The 
sensitivity of a transit survey depends weakly on Rp until 
a critical radius is reached when the S/N of the transit 
falls rapidly. The sensitivity of the survey for 1.0 Rj is 
near this threshold, hence the large scatter in the detec¬ 
tion probability. 

With the detection probabilities for all stars in the sur¬ 
vey for the assumed semimajor axis distribution, we can 
calculate the expected number of detections scaled by the 
fraction of cluster members with planets. Thus, from the 
Poisson distribution, a null result is inconsistent at the 
^95% level when iVdet 3. This allows us to solve for 
the 95% confidence upper limit on the fraction of cluster 
members with planets using Eq. |3 This gives, 

w* 

/*<3.0/^Pdet,i (95%c.l.). (6) 

i^l 

Figure m shows the 95% confidence upper limit on 
the fraction of stars with planets in NGC 1245 for sev¬ 
eral ranges of orbital period. The solid and dashed lines 
give results for 1.5 and 1.0 Rj companions, respectively. 
For 1.5 Rj companions we limit the fraction of cluster 
members with companions to <1.5%, <6.4%, and <52% 
for EHJ, VHJ, and HJ companions, respectively. For 1.0 
Rj companions, we find <2.3% and <15% have EHJ and 
VHJ companions, respectively. 

The detection probability decreases rapidly with or¬ 
bital period beyond 1.0 day. As a result, the survey does 
not reach the sensitivity needed to place an interesting 
upper limit on 1.0 Rj companions beyond P > 3.0 day. 

We further divide the VHJ period range and show 
upper limits for the period ranges l.Oj P/day i2.0 and 
2.0< P/day <3.0, which we denote as P 12 and P 23 - For 
1.5 Rj companions we limit /* to <5.2% and <11% for 
P 12 and P 23 , respectively. For 1.0 Rj companions we 
limit /* to <19% and <47% for P 12 and P 23 , respec¬ 



Fig. 12.— Upper limit (95% Confidence) on the fraction of stars 
in the cluster with companions for several companion radii as label 
along the top. The result for a 1.0 Rj companion is based on the 
entire sample, whereas the results for the other companion radii 
are based on a subsample of TV* = 100 stars. The shaded regions 
denote orbital periods removed by the selection criteria in order 
to eliminate false-positive transit detections that occur around the 
diurnal period and 0.5 day alias. 


tively. We also divide the HJ period range and limit /* 
for 1.5 Rj companions in the 3.0< P/day <6.0 to <36%. 

7.2. Results for Other Companion Radii 

Due to computing limitations we calculate detection 
probabilities for the entire cluster sample only for 1.5 and 
1.0 Rj companions. In 3HI we show that an upper limit 
determination using a subsample of the stars with size 
A^*.SS = 100 approximates the results based on the en¬ 
tire stellar sample. Thus, we calculate upper limits for a 
variety of companion radii using V* sg = 100 randomly 
chosen stars in the sample. Instead of showing upper 
limit results over a range of orbital periods, we derive 
upper limits at fixed period by replacing the semimajor 
axis distribution with a d{a — Oo) function in Equation^ 
To obtain results at fixed period, each star has a different 
Oo that depends on the stellar mass. Figure IT^shows the 
upper limit on the fraction of stars with planets in the 
survey as a function of orbital period. The lines show re¬ 
sults for various values of the companion radius in terms 
of Rj as indicated by the label next to each line along 
the top of the figure. The shaded regions denote orbital 
periods removed by the selection criteria in order to elim¬ 
inate false-positive transit detections that occur around 
the diurnal period and 0.5 day alias. At smaller compan¬ 
ion radii, the transit S/Ncx Rp drops quickly. Toward 
larger companion radii the S/N of the transit saturates 
and the observational window function increasingly dom¬ 
inates the survey effectiveness. The survey cannot detect 
companions with Rp > d.5Rj as the transit/eclipse be¬ 
comes too deep given the removal of measurements that 
deviate by more than 0.5 mag from the mean light-curve 
level. 

8. ERROR IN THE UPPER LIMIT 
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In this section we discuss several sources of error 
present when determining an upper limit on the fraction 
of stars with planets. 


8.1. Error When Using a Subsample 

Computing power limitations discourage calculating 
detection probabilities over the entire cluster sample. 
Thus, we first characterize the error associated with de¬ 
termining an upper limit using only a subset of the entire 
cluster sample. Starting with Equation El we derive an 
error estimate when using a subsample by the following 
means. Replacing the summation over Pi^det with the 
arithmetic mean, (Pjet), Equation El becomes 

/* = 3.0/(iV*(Pdet)). (7) 


By propagation of errors, the error in the upper limit is 
given by 


3.0 <^{p) 

^ ^(Pdet)2’ 


( 8 ) 


where cr^p) is the error in the mean detection probability. 
The error in the mean detection probability scales as 
(T^P^ = crp/yCV*r^, where up is the intrinsic standard 
deviation of the distribution of Pi^det values and A^*,ss is 
the size of the subsample. 

We empirically test this error estimate by calculating 
the upper limit with subsamples of increasing size. The 
small points in Figure 1131 show the upper limit on the 
fraction of stars with planets as a function of the sub¬ 
sample size. The upper limit calculation assumes an even 
logarithmic distribution of semimajor axis for compan¬ 
ions with 1.0 < P < 3.0 day for 1.5 and 1.0 Rj radius 
companions, top and bottom panels, respectively. Neigh¬ 
boring columns of upper limits differ by a factor of 2 in 
the subsample size. We randomly draw stars from the full 
sample without replacement, making each upper limit at 
fixed sample size independent of the others. The dashed 
line represents the upper limit based on the full cluster 
sample. 

The distribution of upper limits around the actual 
value possesses a significant tail toward higher values. 
This tail results from the significant number of stars with 
Ttot = 0.0. At fixed sample size, the large square point 
represents the mean upper limit. Using subsamples sizes 
of iV*^ss ^ 20 tends to systematically overestimates the 
true upper limit. The open star symbol represents the 
1 — (T standard deviation of the distribution at fixed sam¬ 
ple size. The solid line shows the error estimate from 
Equation El Despite the non-Gaussian nature of the 
underlying distribution, the error estimate in the upper 
limit roughly corresponds with its empirical determina¬ 
tion especially toward increasing W,ss where the sys¬ 
tematic effects become negligible. From Figure 113 we 
conclude that adopting IV* ss ^ 100 provides adequate 
control of the random and systematic errors in calculat¬ 
ing an upper limit, without becoming numerically pro¬ 
hibitive. This verifies the procedure for estimating the 
upper limit for a variety of companion radii in a 


8.2. Error in Determining Sample Size 

Up to this point, we have mainly addressed sources 
of error directly associated with determining P^. How¬ 
ever, the upper limit error budget contains an addi¬ 
tional source of error from uncertainties in determining 
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Fig. 13. — Estimates for the upper limit (95% Confidence) on 
the fraction of stars in the cluster as a function of the sample size 
employed in making the estimate are shown as small points. We 
have assumed an even logarithmic distribution in periods between 
1.0 < P/day < 3.0 orbital period. The dash line shows the upper 
limit based on the entire sample. The average upper limit at fixed 
sample size is given by square points. The sequence of open stars 
gives the standard deviation in the distribution of upper limits at 
fixed sample size. The solid line shows the error model estimate 
for the standard deviation in the upper limit. The top panel gives 
results for a 1.5 Rj companion, and the bottom panel gives results 
for a 1.0 Rj companion. 


^mem- Thls additional source of error directly relates to 
the accuracy in determining the number of single main- 
sequence stars in the survey. 

We characterize this error as follows. At fixed orbital 
period, (Pdet) = {PmemPePr)- Given that Pmem is nearly 
independent of the other terms, the previous average is 
separable, such that (Pdet) = {Pmem) {PuPt) ■ This sep¬ 
aration changes the derived upper limit by a negligible 
0.3% relative error. The separation allows us to rewrite 
Equation [3 as 

/<,95 = 3.0/(iV*,eff(P.PT)), (9) 

where fV*^eff = A^*(7mem) is the effective number of clus¬ 
ter members in the sample after taking into account back¬ 
ground contamination. Thus, V*^eff carries equal weight 
with (P^Pt) in the upper-limit error budget. 

The ability to determine iV*^eff accurately provides an 
advantage for transit surveys toward a rich stellar clus¬ 
ter rather than toward a random galactic field. Even 
though methods based on the cluster CMD statistically 
determine cluster membership, they concentrate on a 
narrow main-sequence region to search for planets where 
the cluster counts significantly outweigh the background 
contamination counts. By concentrating on the main se¬ 
quence of a cluster, this survey has only ~ 68% contami¬ 
nation by background stars. In contrast, random galaxy 
fields contain > 90% contami nation by subgiant an d 
giant stars for V<11 surveys llGould fc Morgan! l2003|l . 
Overall, iV*_eff has an 8% error, which propagates to a 
relative error of 8% in the upper limit. The error in 
iV*_eff comes from subtracting the star counts observed 
within a 12.7' radius of the cluster center by the control 
field star counts outside this radius. The error is larger 
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than the Poisson error of = 870 since the control 

field star count is scaled to match the larger cluster field 
area. 


8.3. Error Due to Blends and Binaries 

The final source of error we address results from stellar 
blends due to physical binaries or chance, line-of-sight 
associations. The additional light from an unresolved 
blend dilutes a transit signal from one component of the 
blend. Thus, we overestimate the ability to detect a tran¬ 
sit around blends. However, a compensatory effect arises 
since the extra light from a blend results in an overesti¬ 
mate in the stellar mass and radius, which in turn results 
in modeling a shallower transit. Modeling such details is 
not possible without knowing the binary nature for each 
object, but we can estimate the number of stars affected 
by assuming binary star statistics as measured in the 
field. Due to low stellar crowding, we estimate chance 
blends have a negligible effect in comparison t o physi¬ 
cally associated binaries (IKiss Reddinsll2nnl . Find¬ 
ing charts in Figure El demonstrate the stellar crowding 
conditions of the survey. 

The latest Corav el radial velocity survey dedicated to 
F7-K field dwarfs llHalbwachs et all 1200411 and the vi- 
sual binary and common proper motion pairs survey of 
lEggenberger et all ll2004ll provide the basis for the bi¬ 
nary star estimates. Overall they find a binary frequency 
of 56% for systems with log(P/day) < 6.31. However, 
due to the strong dependence of luminosity on the stel¬ 
lar mass only systems with mass ratio, q > 0.6, signifi¬ 
cantly contribute light to dilute the transit signal. For 
lower mass ratios the lower mass component contributes 
< 20% of the total system flux. When taking binaries 
across the entire range of orbital periods the mass-ratio 
distribut ion peaks near q ^ 0.2 and slowly drops toward 
higher q llDuauennov ^ Ma.vorll 90lTl . From Figure 10 in 


iDiiauennov fc MavoH lIlQQlll . only ^ 20% of their binary 
systems have q > 0.6. Thus, if the binary statistics for 
the cluster matches the field dwarfs, transit dilution oc¬ 
curs for ^ 11% of the stellar sample. The radial velocity 
survey for binaries in the Pleiades and Praesepe clusters 
reveals consis tency with the frequen cy of binaries in the 

In principle, the data from this suryey can also answer 
whether the binary statistics of the cluster matches the 
field dwarfs. Howeyer, the statistical methods and se¬ 
lection criteria described in this study do not optimally 
detect interacting and eclipsing binaries. Additionally, 
in order to reach planetary companion sensitiyities, we 
remoye light-curye deyiations beyond 0.5 mag as dis¬ 
crepant, which remoyes the deep eclipses. 


8.4. Overall Error 

The errors inyolyed with determining the number of 
cluster members dominates the error budget in deter¬ 
mining the upper limit. Howeyer, as discussed in [El 
this is only true if one quantifies and corrects for the 
systematic overestimate in detection probability due to 
a reduction in the transit signal from the procedures of 
generating and correcting the light curve. For instance, 
at the median stellar brightness for this survey, the de¬ 
tection probability is overestimated by >15% for orbital 
periods >4.0 day and >1.0 day for 1.5 and 1.0 Rj com¬ 
panions, respectively, without correction. Since we char¬ 


acterize this systematic effect, the error in determining 
the number of cluster members dominates the error bud¬ 
get. 

Additionally, the potential for a large contamination of 
binaries diluting the transit signal necessitates an asym¬ 
metrical error bar. We roughly quantify the error esti¬ 
mate resulting from binary contamination from the field 
dwarf binary statistics. From the arguments in the pre¬ 
vious section, we adopt ll%asal — cr systematic frac¬ 
tional error due to binary star contamination. Overall, 
combining this systematic error with the 7% fractional 
error in determining the cluster membership, upper lim¬ 
its derived from the full stellar sample contain a 
fractional error. 


9. DISCUSSION 

Along with this work, several other transit sur¬ 
veys have quantified their detection probability from 
actual observations in an attempt to constrain the 
fraction of stars with planets or quantify the consis¬ 
tency with t he solar neighborhoo d radial velocity planet 
disco veries llCilliland et all l200flt iWeldrake et alJ 120051: 

e ieiska et al .1 120051: |Hidas et alJ 1200, 5t iHood et alJ 

. Unfortunately, a direct comparison of upper lim¬ 
its from this work with these other transit surveys cannot 
be made. Until this study, none of the previous studies 
have quantified the random or systematic errors present 
in their techniques in sufficient detail to warrant a com¬ 
parison. Additionally, previous studies do not have quan¬ 
tifiable selection criteria that completely eliminate false¬ 
positive transit detections due to systematic errors in 
the light curve, a necessary component of an automated 
Monte Carlo calculation. 

9.1. Initial Expectations vs. Actual Results 

In the meantime, we can discuss why the initial esti¬ 
mate of finding two planets assuming 1% of stars have Rj 
companions eve nly distributed log arithmically between 
0.03 to 0.3 AU llBurke et alJl200,‘^ compares to the re¬ 
sults from this study, which indicate that we expected 
to detect 0.1 planets. The initial estimates for the de¬ 
tection rate are based on the theoretical framework of 
ICaiidil ll2000ll . Civen a photometric noise model, obser¬ 
vational window, and S/N of the transit selection crite¬ 
ria, the theoretical framework yields an estimate of the 
survey detection probability. This theoretical detection 
probability coupled with a luminosity function for the 
cluster determines the expected number of detections. 
As we show next, the initial estimates did not account 
for the light curve noise floor or detector saturation, and 
contain optimistic estimates for the sky background and 
luminosity function. In addition, the initial estimates 
could not have accounted for the 50% reduction in sig¬ 
nal for the majority of th e lig ht curves due to the de¬ 
tector er ror discussed in [lliHI Finally, as discussed in 
detail bv iPepoer fc Caudil ll2005ll and demonstrated ex¬ 
plicitly here, the detection probability is very sensitive 
to the precise error properties near the critical threshold 
of detection, which for this survey is just reached for Rj 
companions. 

The top panels of Figure El compare the detection 
probability of the Monte Carlo calculation of this study 
to the initial theoretical estimate. The small points repli¬ 
cate the Monte Carlo results from the top panels of Fig- 
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Fig. 14. — Top-. Probability for transit detection as a function of 
the apparent 7-band magnitude assuming an even logarithmic dis¬ 
tribution in semimajor axis from 1.0< P <3.0 day and Pmem = 1-0 
using the Monte Carlo calculation of this study(sma/^ points). The 
binned average of the Monte Carlo results is denoted by open stars. 
The dash line shows the expected probability for transit detection 
based on a theoretical calculation prior to this survey. The dot dash 
line shows the theoretical probability for transit detection assum¬ 
ing a photometric noise model appropriate for the survey. The solid 
line shows the theoretical probability for transit detection with an 
accurate photometric noise model for the survey and including the 
effects of limb darkening. The left panel shows 1.5 Rj companion 
results. The panel shows 1.0 Rj companion results. Bottom: 

Shows the theoretical probability for transit detection allowing each 
star of the survey to have its empirically determined photometric 
noise and including the effects of limb darkening {small points). 
The open stars are reproduced from the top panels. 


ure cni while the dashed line shows the detection proba¬ 
bility based on the initial theoretical expectations. The 
initial theoretical expectations clearly overestimate the 
detection probability. The bright end continues to rise 
due to ignoring the effects of detector saturation and the 
photometric noise floor. The faint end does not cutoff 
due to an underestimated sky brightness. The initial 
estimate of the sky brightness, 19.5 mag arcsec“^, com¬ 
pares optimistically to the range of sky brightnesses en¬ 
countered during the actual observations. The sky varied 
between 17.5 and 19.0 mag arcsec“^ over the course of 
the observations. The full lunar phase took place near 
the middle of the observation, and the Moon came within 
40° of the cluster when nearly full. 

The initial estimate for the cluster luminosity func¬ 
tion simply selected cluster members via tracing by eye 
lines that bracket the main sequence in the CMD. This 
crude technique led to an estimated 3200 cluster mem¬ 
bers down to I ^20. A careful accounting of the field 
star contamination results in only ~870 cluster mem¬ 
bers in the survey. The luminosity function overestimate 
and the expected sensitivity to transits around the bright 


and faint cluster members leads to a factor of 4-5 overes¬ 
timate in the number of cluster members in the survey. 
Additionally, the factor of 4-5 overestimate of the initial 
detection probability when compared to binned average 
detection probability for the Monte Carlo results (open 
stars in Figure nui . easily accounts for the factor of 20 
difference in the overall number of expected detections 
(for R = Rj). 

9.2. Improving Theoretical Expectations 

Clearly, accurate and realistic transit detection statis¬ 
tics requires more detailed analysis than these early es¬ 
timates an d more careful theoreti cal work has already 
been done l|PeDDer fc Caudil 1200511 . In the case of an 
open cluster, delineating cluster membership by tracing 
the main sequence in the CMD overestimates the num¬ 
ber of cluster members. A careful subtraction of the field 
contamination is necessary in order to extract an accu¬ 
rate cluster-member count. 

A photometric noise model that accurately reflects the 
quality of observations is the next step in correctly calcu¬ 
lating a theoretical detection probability. From Figure^ 
we estimate the actual photometric noise present in the 
data. This includes the proper sky measurement and 
systematic floor in the photometric precision. With a 
noise model similar to the lower solid line in Figure ^ 
we recalculate the theoretical detection probability. The 
dot dash line in Figure m shows the resulting detection 
probability still overestimates the Monte Carlo results. 
However, it does agree with the faint-end cutoff of the 
Monte Carlo calculation. We impose the bright-end cut¬ 
off due to saturation effects at the same magnitude as the 
observed increase in light curve rms as shown in Figure^ 

For these results we i nclude a n add itional effect not 
taken into account by IG^ We m ultiply 

the transit S/N selection criteria. Equation 5 ofiGaudil 
(j2QQ3): by -\/max(iVobs, 1-7), where Vobs is the typical 
number of transits detected throughout the observing 
run. The iVobs = 1-7 floor in this factor corresponds 
to the requirement of observing the transit twice multi¬ 
plied by the observing efficiency. For simplicity, we take 
Noha = Ntot/P X 0.2, where Wot = 16, the length of the 
observing run in days, and the factor of 0.2 accounts for 
the actual observational coverage encountered during the 
run. 

Given that the theoretical calculation still overesti¬ 
mates the Monte Carlo results, to increase the realism 
of the theoretical detection probability, we include a lin¬ 
ear limb-darkening law, which effectively weakens the 
tr ansit depth. We solve for the factor G, Equation 6 
ofiGaudil i|2000l) . assuming a linear limb-darkening pa¬ 
rameter, fj, = 0.6, for all stars. The inclusion of limb 
darkening significantly impacts the theoretical detection 
probability as the solid line in Figure m demonstrates. 
Although the theoretical detection probability still over¬ 
estimates the upper envelope of results from the Monte 
Carlo calculation, the level of agreement, after including 
an accurate photometric noise model and limb darkening, 
shows significant improvement over the initial estimates. 

Despite the improved agreement, the Monte Carlo de¬ 
tection probability calculation shows significant scatter 
at fixed magnitude. The theoretical probability treats 
all stars at fixed magnitude as having the same noise 
properties. With the theoretical detection probability we 
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can address whether the scatter in detection probability 
at fixed magnitudes results from the observed scatter in 
noise properties at fixed magnitude as shown in Figure ^ 
Thus, we calculate a theoretical detection probability for 
each star individually using the measured rms in the light 
curve for each star to determine the the oretical transi t 
S/N selection criteria using Equation 5 of iGaudil ll200f]ll . 
The small points in the bottom panels of Figure [HI show 
the resulting theoretical detection probability. 

Some of the scatter in detection probability results 
from the scatter in noise properties as a function of 
magnitude. The heavy star points represent the aver¬ 
age Monte Carlo detection probability in 0.25 magnitude 
bins. In the case of the 1.5 Rj companions, the signal is 
large in comparison to the photometric noise. The left 
panels of Figure El demonstrate the theoretical detec¬ 
tion probability overestimates the Monte Carlo detection 
probability by only 20%. However, the closer the transit 
signal approaches the systematic and rms noise, the the¬ 
oretical detection probability strongly overestimates the 
actual detection probability. In the case of 1.0 Rj com¬ 
panions (right panels of Figure m, the theoretical cal¬ 
culation overestimates the Monte Carlo results by 80%. 
Thus, we urge caution when relying on a theoretical de¬ 
tection probability when the survey is near the critical 
threshold for transit detection. Such is the case for 1.0 
Rj companions in this survey. 

9.3. Planning Future Surveys 

Even though the theoretical calculation overestimates 
the absolute detection probability by a factor of <2, tests 
on a small sample of stars with the Monte Carlo calcu¬ 
lation reveal it provides a much higher relative accuracy. 
Thus, the computationally efficient theoretical calcula¬ 
tion allows us to examine the relative change in the de¬ 
tection probability for a given change in survey parame¬ 
ters. For planning future surveys it is essential to decide 
between increasing the number of stars by observing an¬ 
other cluster or improving the detection probability by 
increasing the length of observations on a single cluster. 
As shown in the upper limit scales linearly with the 
sample size, thus keeping everything else constant, in¬ 
creasing the sample size by a factor of 2 improves the 
upper limit by a factor of 2. 

Using the theoretical detection probability framework, 
we can quantify the improvement in sensitivity for a sur¬ 
vey twice as long. We assume a survey twice as long 
consists of an observing window identical to the current 
survey for the first half and repeats the observing win¬ 
dow of the current survey for the latter half. The upper 
limit improves only by a factor of 1.3 for a logarithmic 
distribution of VHJ planets. However, the upper limits 
for HJs with 3.0 to 9.0 day orbital periods decrease by a 
factor of 2.6. Thus, not only is it more efficie nt to observe 
this cl uster twice as long, but the analysis of iCaudi et al.l 
<I2005^ reveals a 5-10 times larger HJ population than the 
the VHJ population. This strongly suggests transit sur¬ 
veys with a single observing site require month long runs 
for maximum efficiency in detecting HJ companions. 

Figure m reveals little improvement in the detection 
probability occurs for increasing the photometric pre¬ 
cision, at least for 1.5 Rj companions. To first order, 
the photometric precision determines the faint-end cut¬ 
off in the detection probability. Thus, a lower sky back¬ 


ground or improved photometric precision predominately 
effects the number of stars in the survey rather than 
the detection probability. However, improving the pho¬ 
tometric precision does lead to increasing the sensitiv¬ 
ity for smaller radius companions. In the case of 1.0 
Rj companions, the rms in the light curve typically is 
<1.8 times lower than the transit signal. As shown in 
the previous section, the theoretical detection probabil¬ 
ity breaks down for such low precision. In the case of 
1.5 Rj companions, the rms in the light curve typically 
is <4 times lower than the transit signal. Thus, for the 
I.O Rj results to reach the same sensitivity as the 1.5 
Rj results, improvement in the light curve rms is nec¬ 
essary until the transit S/N is above a critical threshold 
whe n the detection probab ility is weakly dependent on 
Rp llPepper fc Caudiil2005|l. 

According to a recent review of radial velocity detected 
planets, 1.2 ± 0.3% of solar neig hborhood stars have HJ 
companions llMa.rc.v et 3,1.112(1(1^ . This survey of NCC 
1245 reached an upper limit of 52% of the stars having 
1.5 Rj HJ companions. As mentioned previously, a sur¬ 
vey lasting twice as long can reduce this upper limit to 
21%. Reaching similar sensitivity as the radial velocity 
results requires observing additional clusters in order to 
increase the number of stars in the sample. This sur¬ 
vey has ~ 870 cluster members and ~ 740 of them have 
nonzero detection probability for 1.5 Rj VHJ compan¬ 
ions. Hence a total sample size of ^ 7400 dwarf stars ob¬ 
served for a month will be needed to help constrain the 
fraction of stars with planets to a 2% level (comparable 
to radial velocity results). Assuming that the observed 
HJ frequency of ~ 1% remains valid for a variety of stel¬ 
lar environments, we expect to detect one planet every 
5000 dwarf stars observed for a month. Results for 1.0 
Rj companions without substantial improvement in the 
photometric precision likely will require a small factor 
larger sample size. 

10. CONCLUSION 

In this study we complete the analysis of a 19-night 
search for transiting extrasolar planets orbiting mem¬ 
bers of the open cluster NCC 1245. An automated tran¬ 
sit search algorithm with quantitative selection criteria 
finds six transit candidates; none are bona fide planetary 
transits. Thus, this work also details the procedure for 
analyzing the null-result transit search in order to deter¬ 
mine an upper limit on the fraction of stars in the cluster 
harboring close-in Rj companions. In addition, we out¬ 
line a new differential photometry technique that reduces 
the level of systematic errors in the light curve. 

A reliable upper limit requires quantifiable transit se¬ 
lection criteria that do not rely on visual, qualitative 
judgments of the significance of a transit. Thus, we 
develop completely quantitative selection criteria that 
enable us to calculate the detection probability of the 
survey via Monte Carlo techniques. We inject realistic 
limb-darkened transits in the light curves and attempt 
their recovery. For each star we inject 100,000 transits 
at a variety of semimajor axes, orbital inclination angles, 
and transit phases, to fully map the detection probability 
for 2700 light curves consistent with cluster membership 
based on their position in the CMD. After characterizing 
the field contamination, we conclude the sample contains 
~870 cluster members. 
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When calculating a 95% confidence upper limit on the 
fraction of stars with planets, we assume companions 
have an even logarithmic distribution in semimajor axis 
over several ranges of orbit al period. W e adop t the pe¬ 
riod ranges as outlined bv iGaudi et al.l ll2005ll . for HJ 
and VHJ companions, and an as of yet undetected pop¬ 
ulation with P<1.0 day, which we denote as Extremely 
Hot Jupiters (EHJ). For NGC 1245, we limit the fraction 
of cluster members with 1.0 Rj companions to <3.2% 
and <24% for EHJ and VHJ companions, respectively. 
We do not reach the sensitivity to place any meaningful 
constraints on 1.0 Rj HJ companions. For 1.5 Rj com¬ 
panions we limit the fraction of cluster members with 
companions to <1.5%, <6.4%, and <52% for EHJ, VHJ, 
and HJ companions, respectively. 

We also fully characterize the errors associated with 
calculating the upper limit. We find the overall error 
budget separates into two equal contributions from er¬ 
ror in the total number of single dwarf cluster members 
in the sample and the error in the detection probability. 
After correcting the detection probability for systematic 
overestimates that become increasingly important for de¬ 
tecting transits toward longer orbital periods (see , we 
conclude that random and systematic errors in determin¬ 
ing the number single dwarf stars in the sample dominate 
the error budget. m details the error analysis, and over¬ 
all, we assign a ^ 7 ^ fractional error in the upper limits. 

In planning future transit surveys, we demonstrate 
that observing NGC 1245 for twice as long will reduce 
the upper limits for the important HJ period range more 
efficiently than observing an additional cluster of similar 
richness as NGC 1245 for the same length of time as this 
data set. To reach a ^ 2% upper limit on the fraction of 
stars with 1.5 Rj HJ companion s, where radial velo city 
surveys currently measure 1.3% llMarcv et al.ll200,'i^ . we 
conclude a total sample size of ^ 7400 dwarf stars ob¬ 
served for a month will be needed. If 1% of stars have 1.5 
Rj HJ extrasolar planets, we expect to detect one planet 
every 5000 dwarf stars observed for a month. Results for 
1.0 Rj companions without substantial improvement in 
the photometric precision likely will require a small fac¬ 
tor larger sample size. 
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Fig. 15.— The lines show the difference in between a 

boxcar fit to a planetary transit across a limb-darkened star and the 
exact model fit, normalized by Axq the difference in x^ between 
the exact model fit and a constant fiux fix to the light curve. Each 
band is for a different planet/star radius ratio Rp/Rt, and the 
width of the band shows the variation in Ax^^^/Axq for range of 
linear limb-darkening parameters ui = 0.0 — 0.4. 
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APPENDIX 

In this appendix, we derive the boxca r fitting algorithm used to search for planetary transits. The original developers 
of this algorithm llKovacs et a,l .1 l2fi?)^ . hereafter BLS, study its properties in the presence of Gaussian noise. We 
rederive the BLS method he re in terms of the fami liar language of minimization, demonstrating the equivalence of 
the two approaches (se e also lAigrain fc Irwinl l2004ll . In the process, we elucidate the meaning of some of the detection 
statistics introduced in lKovac^^^ ll2002ll . Finally, we quantify the accuracy of the boxcar approximation to transit 
light curves and discuss appropriate detection thresholds in the presence of correlated systematic errors as well as pure 
noise. 

We begin by considering a series of flux measurements, Fk, taken at times tk with errors (TF.fe- The equation, 
^ — (Sfe ^k/o'pj^)/{J2k defines the mean error-weighted flux of the measurements. We can write the fractional 

deviations from the mean as ruk = F^/F — 1 and the fractional errors as ak = <JF,kl F. For small deviations, the 
difference from the mean in magnitudes approximates — 1.086mfc. 

A light curve of a star with transiting planet can be written as F{t) = Eo[l — J(t)], where Eq is the unobscured flux 
of the star, and 5{t) is the fractional drop in flux due to the transiting planet as a function of time. This generally 
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depends on five parameters, namely the planet period P, the phase of the planet orbit the time it takes for the 
planet to cross one stellar radius, the ratio of the planet radius to stellar radius p, and the limb-darkening parameter 
of the star. We can write the fractional flux difference as m{t) = F{t)/F — 1 = mo — 5{t) — moS{t), where we have 
defined mo = Fq/ F — 1. We will assume that i5 <C 1 and mo 1, and therefore m{t) ~ mo — S{t). 

For p 1, and no limb darkening, the light curve m{t) phased to the period P reduces to a simple boxcar, 5 = 
when the planet is occulting the star and i5 = 0 otherwise, with five parameters: the unocculted flux mo, the transit 
period P, phase (f), duration At, and depth 6. This simple form allows an analytical solution for two of these parameters. 
For finite p and limb darkening, the expression for 6{t) increases in complexity, but, as we show later, a simple boxcar 
still approximates the variability. 

The likelihood C that the model m(t) with parameters a = (mo, 6, P, At) describes the data ruk is simply, 


— 2\n£ = 
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mfc - m(tfc) ] ^ 

O'fc 


( 1 ) 


A likelihood maximization, or minimization, determines the best-fit parameters. Phasing the measurements mk to 
a given period, we can split into two terms, including the points in transit and points out of transit. 
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where the first sum over index i is over the points in transit, and the second over index j is over the points out of 
transit. The quantity P of BLS corresponds to our X^/X^fe where the sum over index k is over all points. Since 
Sfe is a constant, this verifies the equivalence of their algorithm to a simple x^-minimization. Expanding the 

quadratic terms in equation @ yields, 
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Given the identity m = {J2k = 3, it is clear that the first term is simply the x^ of a constant flux 

fit to the data, which we denote Xo^ the second term is zero. Therefore, the last four terms are the improvement 
in x^ between the constant-flux fit, and a boxcar transit fit, for a given P, (j), At. 

Minimizing the expression for x^ with respect to mo, we find mg = S/cInserting this into the 

expression for x^, and then minimizing with respect to <5, we find the parameters <5, mg that simultaneously minimize 
X^ for a given P, cj), At, 
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We note that the solution for 5 does not impose a particular sign. The best-fit boxcar model may be of a transit or 
anti-transit nature. The improvement in x^ at these best-fit parameters is 
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At fixed P, (j), At, and for pure Gaussian noise, Ax^ is distributed as x^ with one degree of freedom (corresponding to 
the one additional free parameter). 

The above expression gives the x^ improvement for a particular set of P, (j). At. The global solution that maximizes the 
X^ improvement, Axmaxj requires a grid search over the entire P, (p. At parameter space feasible for transit detection. 
In the BLS algorithm, evaluation of Ax^ for a given P, p, At amounts to simple error-weighted binning, making the 
process extremely efficient and fast. Due to the large number of light curves with injected transits that must be searched 
in order to determine the detection probability accurately (see this study requires the excellent numerical efficiency 
provided by the BLS method. Comparison with BLS reveals that Ay^ax has a close correspondence with their detection 
statistic SR, such that Ax^ax/Efc ^k^ = —(SR)^. If we assume constant noise, i.e. ak = constant = a, and define Nt 

to be the number of points in transit, then the effective signal-to-noise ratio is S/N = (jAxmaxI)^^^ = 

Although, at fixed P, (j>. At, analytic expressions can provide the significance of a given value of x^ in the presence 
of pure Gaussian noise, in the present case folding data over many trial periods and phases hampers determination 
of the significance of the global Ax^ax- Thus, assessing the probability of achieving the observed outcome by chance 
involves taking into account the effective number of independent trials. Several papers (e.g., BLS. IJenkins et ffi]l2002ll 
propose methods to address these issues. However, as we have shown, in the presence of correlated systematic errors, 
formulations with a Gaussian-noise basis underestimate the detection thresholds needed to avoid falsely triggering on 
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the systematic errors. Generally, one must use the properties of the data themselves to set the appropriate detection 
threshold. 

In general, a matched filter optimally detects a signal of known form in noisy data. However, the efficiency of boxcar 
fitting makes it highly advantageous, and as we show next nearly optimal for transit detection. This is because, for 
small planets and modest limb darkening, simple boxcars match well the planetary transit signal. We demonstrate this 
in FigureEl where we show = Xbc ~ Xexact^ the difference in between a boxcar fit to a planetary transit across 
a limb-darkened star and the exact model fit, assuming uniform errors. In order to make the result independent of the 
error properties of the light curve, we plot this normalized to Ayp = xionstant ~ Xexact^ ftie difference in y^ between a 
constant flux fit to the light curve and the exact model fit . In other words, since (Ay^)^/^ is the total signal-to-noise 
ratio of the transit detection, the ratio Ay^^/Ayp is (the square of) the fractional degradation of the signal-to-noise 
ratio that arises from fitting the approximate boxcar form to the true transit shape. We show Ay^^/Ayp versus the 
impact parameter of the transit for four different values of p = 0.01, 0.05, 0.1 and 0.2, as well as linear limb-darkening 
coefficients in the range ui = 0 — 0.4. The fractional difference in Ayg is less than ~ 5% for impact parameters < 0.8 
and p < 0.1, and it is always < 10% for the range of parameters relevant here. However, if the number of detections is 
a strong function of y^ near the detection threshold, then even a ~ 5% change in y^ can have a significant effect on 
the inferred detection probability. Therefore, it is important to inject realistic transits when determining the recovery 
rate and detection probability, to account for the imperfect match of the boxcar filter. 
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Table 1. MDM 2.4m Observations 


Date (2001) 

#Exps 

FWHM (arcsec) 

Comments 

Oct. 24 

75 

1.2 

clear - 1st quarter moon 

25 

73 

1.4 

partly cloudy 

26 

67 

1.4 

partly cloudy 

27 

22 

1.6 

overcast 

28 

96 

1.4 

cirrus 

29 

86 

1.3 

cirrus 

31 

32 

1.4 

partly cloudy 

Nov. 1 

32 

1.4 

clear humid - full moon 

2 

80 

1.5 

clear - moon closest approach 

6 

44 

1.4 

clear humid 

7 

92 

1.3 

partly cloudy - 3rd quarter moon 

8 

57 

1.4 

cirrus 

10 

81 

1.8 

cirrus 

11 

99 

1.3 

clear 



Table 2. Transit candidate data 


to 


ID 

RA(2000.0) 

Dec(2000.0) 

V (mag) 

B — V(mag) 

V — I (mag) 

Xmem 

P(d) 

A/(mag) 

r(h) 

0 

Ax^/Ax?. 

Ax^ 

/ 

M(Mq) 

Log((R/R©) 

Teff(K) 

30207 

3:15:40.0 

+47:21:18 

18.1 

1.17 

1.31 

0.137 

4.614 

0.030 

1.66 

0.72 

5.65 

584 

0.57 




20513 

3:15:04.6 

+47:15:09 

18.6 

1.10 

1.27 

0.028 

1.637 

0.018 

4.71 

0.95 

3.85 

840 

0.49 

0.91 

-0.095 

5400 

20065 

3:15:03.8 

+47:14:33 

16.1 

1.02 

1.25 

0.418 

3.026 

0.115 

4.18 

0.24 

3.82 

88086 

0.59 




20398 

3:14:49.5 

+47:16:03 

18.4 

1.29 

2.00 

0.863 

0.349 

0.145 

0.90 

0.72 

3.76 

66056 

0.29 




20274 

3:14:35.9 

+47:19:29 

19.3 

1.67 

3.37 

4.390 

0.302 

0.050 

0.98 

0.48 

3.23 

14063 

0.21 




70718 

3:13:56.2 

+47:06:59 

21.1 

1.28 

1.79 

0.017 

0.640 

0.032 

3.07 

0.16 

2.86 

312 

0.22 

0.63 

-0.253 

4300 
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